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SIMPLIFIED PROCEDURES AND CHARTS FOR THE RAPID ESTIMATION 
OF BENDING FREQUENCIES OF ROTATING BEAMS^ 

By Robert T. Yntema 


SUMMARY 


A Rayleigh energy approach utilizing the bending mode of the nonro- 
tating beam in the determination of the bending freq^uency of the rotating 
beam is evaluated and is found to give good practical res^llts for heli- 
copter blades . 

Charts are presented for the rapid estimation of the first three 
bending frequencies for rotating and nonrotating cantilever and hinged 
beams with variable mass and stiffness distributions^ as well as with 
root offsets from the axis of rotation. Some attention is also given to 
the case of rotating beams with a tip mass . 

A more exact mode -expansion method used in evaluating the Rayleigh 
approach is also described. Numerous mode shapes and derivatives obtained 
in conjvinction with the frequency calculations are presented in tabular 
form. 


INTRODUCTION 


Designers of helicopter rotor blades generally agree that accurate 
means are needed for estimating the natural bending frequencies of the 
rotating blades in order to obtain a blade design which is as free as 
possible from resonant or near -resonant excitation by the periodic 
loading on the rotor. Although numerous methods are available for deter- 
mining the bending frequencies of rotating blades (see, for example, 
refs. 1 to 1^+), designers have expressed the need for a simplified, yet 
reasonably accurate, procedure for their determination, preferably in 
the form of a set of charts. With this need in mind, an investigation 
was undertaken which had a twofold purpose: (a) an evaluation to show 

whether a Rayleigh energy approach utilizing the mode shape of the non- 
rotating beam may be employed to obtain close approximations for the 
natural bending frequencies of the rotating beam and (b) a set of charts 

^An amplified and extended version of NACA PIM L5^G02, 195^* 
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which would permit the rapid estimation of the first three bending fre- 
quencies of both nonrotating and rotating hinged and cantilever blades. 

The main purpose of this report is to present this evaliiation and the 
frequency charts. 

The Rayleigh energy approach was evaluated with respect to such 
items as various rotational speeds, higher modes, flapping-hinge or root 
offset, variable blade mass and stiffness distributions, and a large con- 
centrated tip mass. The evaluation was made by comparing frequency 
results obtained by the Rayleigh method with results obtained by a more 
accurate mode -expansion method. The details of the mode -expansion method 
are given in appendix A. 

The charts for frequency estimation were obtained by considering 
various families of beams with selected mass and stiffness distributions 
and were derived for both hinged and cantilever beams. The frequencies of 
both nonrotating and rotating cases may be estimated for (a) beams with 
and without offset which have mass and stiffness distributions which can 
be approximated by linear relations and (b) beams with uniform mass and 
stiffness distributions plus a concentrated mass at the tip. 

If the bending frequencies of the nonrotating beams are known, a 
third set of charts permits the estimation of the bending frequencies of 
rotating beams with approximately linear stiffness distribution and arbi- 
trary mass distribution. 

As an adjunct to the Rayleigh approach utilizing the nonrotating- 
beam mode shapes, a method is presented in appendix B which permits a 
fairly accurate determination of the first bending mode and frequency of 
a rotating or nonrotating hinged beam with a tip mass from a knowledge 
of the first bending mode of the nonrotating beam without a tip mass . 

The report also presents bending -mode results, obtained in conjunc- 
tion with the frequency determination, which show the effect of the parajn- 
eters on mode shape. Many of these mode shapes are tabulated in normal- 
ized form together with their first and second derivatives, or as mode 
coefficients (coefficients of an expansion in terms of uniform-beam modes). 
These results can be used in connection with the modified approach of 
appendix B or in other analyses . 

In order to facilitate the further application of the mode -expans ion 
method to the accurate determination of modes and frequencies of rotating 
beams with linear mass and stiffness distributions, concentrated tip mass, 
and offset different from those considered herein, certain integrals which 
have been evaluated are also presented in tabular foirni. 
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SMBOLS 

An 

<1 

mode coefficients for the nth rotating-beam mode 
(see eq. (A3)) 

a 

1 ^ 

nonrotating-beam bending frequency coefficient, — — 

V EIq 

b 

/nuL^ 

rotating -beam bending frequency coefficient, — 

V ^^0 

c 

beam-stiffness -distribution constant (see eq. (a12)) 


pendulum- or zeroeth-mode coefficient (see appendix B) 

El(x) 

leng^hwise bending stiffness distribution for beam 

EIo 

bending stiffness of beam at root 

Ei(x) 

nondimens ional bending stiffness distribution for 
beam, EI(x)/eIq 

K 

Southwell coefficient (see eqs. (4) and (5)) 

Ko 

zero-offset Southwell coefficient 

Kl 

offset-correction factor for Southwell coefficient 

Ko’ 

zero-offset rotating-beam frequency coefficient} found to 
be essentially independent of beam mass distribution 
(see eq. (ll)) 

k 

beam-mass -distribution constant (see eq. (A12)) 

L 

beam length, meas\ired from point of root fixity to tip 

m(x) 

lengthwise mass distribution for beam (mass per unit length) 

mo 

mass per unit length of beam at root 

m(X) 

nondimensional mass distribution for beam, m(x)/mo 

nid 

part of beam mass distribution which is continuous (not 
concentrated) 

Mt 

mass concentrated at tip of beam 


k 
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r 


T 



X 

X 

Y 

y 

5(x-L) 

6(x-l) 

e 

e 

0 


0 

fi 

“NR 


nondimensional tip-mass ratio, 

lengthwise distribution of tension force in beam, Ti(x)ii2 

lengthwise-distribution function for tension force 
(see eq. (2) ) 

spanwise coordinate along beam measiured from root 
nondimensional spanwise coordinate, x/l 
bending mode shape of nonrotating beam 
bending mode shape of rotating beam 
Dirac delta function 

Dirac delta function in nondimensional coordinates 
offset of hinge or point of fixity from axis of rotation 
nondimensional offset, ejh 
dummy variables for x and x 

characteristic nimiber for nonrotating uniform beam with 
mass at tip; identical to sq\iare root of nonrotating- 
beam bending frequency coefficient 

bending mode shape for nonrotating uniform beam normalized 
at tip 

rotational speed of beam 

natural bending frequency of rotating beam 
natural bending frequency of nonrotating beam 


Subscripts : 

n integral number designating natural bending mode of beam 

F beam cantilevered or fixed at root 

t tip of beam 

Primes mean differentiation with respect to x or x unless 
indicated otherwise. 
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THE RAYLEIGH APPROACH APPLIED TO A ROTATING BEAM 

Description 


The problem being treated in this report is a rotating beam vibrating 
freely in one of it.s natural bending modes. By equating the kinetic energy 
at zero displacement to the potential energy of both the bending and cen- 
trifugal forces at maximum displacement, the following frequency equation 
for vibration perpendicular to the plane of rotation can easily be derived: 



where n refers to the mode under consideration and 


"^ 2 . = + e)m dq (p) 

Jx 

Equation (l) yields an exact value for the nth bending frequency of a 
beam rotating at any rotational speed G if the nth natural bending mode 
shape of the rotating beam is known for this value of G. Unfortunately, 
the mode shape is usvially just as much of an unknown as the frequency is. 
An estimation of the frequency may be made, however, by making use of the 
well-known Rayleigh principle; that is, a mode shape which is consistent 
with the constraints of the system is assumed and is used to evaliiate the 
energy integrals which, in turn, give an approximate value for the fre- 
quency. In this report the nonrotating -beam mode shape is chosen as the 
approximation for the rotating -beam mode shape, and an evaluation is made 
to show whether the use of such a shape yields close approximations to 
the exact frequencies of the rotating beam. 

If the nth mode shape of the nonrotating beam Yn is substituted 

into equation (l), the first term becomes exactly the square of the 
bending frequency of the nonrotating beam. By denoting the ratio in the 
second term by a Southwell coefficient, the frequency equation takes 

the following simplified form: 




( 5 ) 
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where 


Kn = 



^dx 


■L 

(t) + e)m dT) 



dx 


(^) 


This expression for can be subdivided into two independent parts 

as follows: 


Kn = X 


K-i e 
“^n 


(5) 


where both 
defined as 


Kn and 
follows : 



are independent of the offset 


e 


and are 



( 6 ) 


In the remainder of this report is referred to as the zero-offset 

Southwell coefficient and Ki is referred to as the offset-correction 

“^n 

factor for the Southwell coefficient. 

2 

It is convenient to write in terros of a nonrotating-beam 

frequency coefficient and the mass and stiffness of the beam at the 

root as 
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(7) 


By means of equations (5) and ( 7 ), equation ( 5 ) may be written as 


Equation ( 8 ) with Kq , % ® nondimens ional form serves 

as the basis for the charts for rapid frequency estimation to be pre- 
sented subsequently in this report. These charts pro->/ide values of a^, 

Kn > and Kt which, in conjunction with the mass and stiffness of the 
^n -‘-n 

beam at the root, the length of the beam, the hinge offset, and the rota- 
tional speed, permit rapid estimation of the first three bending fre- 
quencies of rotating or nonrotating beams . 

If the mass distribution of the blade is given by a simple analytic 
function, the integral expression for (eq. ( 2 )) can usually be 

evaluated exactly; for arbitrary mass distribution, however, numerical- 
integration methods such as are given in reference I 5 must be employed. 
Because of the natiare of the numerical. -integration procedirre used in the 
present paper, a slightly different form of the expression for Kn was 

found to be usefilL. This form can be obtained by performing an integra- 
tion by parts on the mmierator of equation (4), and in nondimens ional 
form the result appears as 




(9) 


whence the definitions for Kn and K-, are evident. 

n ->-n 


An additional form of equation ( 3 ) is now presented for use in sub- 
sequent sections of this report. Dividing equation ( 5 ) by yields 
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This form of equation (3) was foimd to be useful in the evaluation of the 
Rayleigh approach. Hereinafter, in this report is referred to 

as the frequency parameter and is referred to as the rotational- 

speed parameter. Also, for subsequent use in this report, a new zero- 
offset rotating-beam frequency coefficient defined as 



where the subscript F indicates that a2^ is the nonrotating -beam 

frequency coefficient for the beam cantilevered at the root. All other 
terms are for the beam with its actual root condition, that is, either 
cantilevered or hinged. 

It is shown subsequently in this report that this new frequency 
coefficient is insensitive to beam mass distribution and should there- 
fore be iisef\il in estimating bending frequencies for families of beams 
with similar stiffness distributions. As is apparent from equation (ll), 
the fundamental frequency of the nonrotating beam treated as a cantilever 
must be known in addition to the bending frequencies of the beam with 
the actual root condition (cantilevered or hinged). 


Evaluation of Rayleigh Approach 

In order to determine the accuracy, usefulness, and possible limita- 
tions of the Rayleigh approach based on nonrotating -beam bending modes, 
the bending frequencies were calculated by this approach for a series of 
rotating beams with systematically varied parameters; the frequencies 
obtained in this manner are compared in this section with the results 
obtained by the more exact mode -expans ion method of appendix A. For the 
cantilever beams, five uniform-cantilever -beam bending modes were used in 
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the expansion; for the hinged beams ^ a pendulum mode was included in 
addition to five hinged-beam bending modes. 

The cases studied by both methods are shown in figure 1. Both 
cantilever and hinged beams are considered for the following cases: 

(a) Uniform beams with 0- and 10-percent root offset 

(b) Beams with mass and stiffness distributions varying linearly 
from the root value to zero at the tip and with 0- and 10-percent root 
offset 

(c) Uniform beams with a mass at the tip. 

The results for all the cases treated were obtained in nondimen- 
sional form and are presented in plots in which the variation of bending 
frequencies with rotational speed as predicted by the exact method of 
appendix A and by the Rayleigh approach may be compared. In each of the 
figures introduced in this section the abscissa is the squared nondimen- 
sional rotational-speed parameter (the squared ratio of rotational speed 
to the first bending frequency of the nonrotating beam) and the ordinate 
is the squared nondime ns ional frequency parameter (the squared ratio of 
the nth bending frequency of the rotating beam to the nth bending fre- 
quency of the nonrotating beam). 

The range of the rotational-speed parameter in each case corresponds 
roughly to that encountered in current helicopters with some latitude for 
new design. Since the first bending frequency of a hinged beam is roughly 
four times the first bending frequency of the same beam fixed at the root, 
widely different scales result for the hinged and cantilever beams. The 
abscissa range also varied with tip mass because the fundamental fre- 
quency of a nonrotating beam decreases with increase in tip mass. For 
the -uniform cantilever beam with a tip mass, this variation is large and 
thus results in a greatly expanded abscissa scale with each increase in 
tip mass. For the uniform hinged beam with a tip mass, the effect of 
tip mass on the nonrotating frequency is relatively small and thus the 
abscissa range was not extended appreciably with each increase in tip 
mass . 


Hinged beams without tip mass .- The variation of bending frequency 

with rotational speed for a uniform* hinged beam is shown in figure 2 for 
offsets of 0 and 10 percent. For this case the Rayleigh approach may be 
seen to be very accurate for all three modes throughout the entire 
rotational-speed range covered. The maximum error is about 5 percent in 

the frequency squared and thus only about li percent in the frequency. 

2 

This maximum error occurs at the highest rotational speed and is roughly 
the same for all three modes. 
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Frec^uency results for the hinged heads with linear mass and stiff- 
ness distributions are shown in figure 5 for offsets of 0 and 10 percent. 
From this figure it is evident that the results obtained by the Rayleigh 
method for this case are very accurate, even for the highest rotational 
speeds shown. 

A comparison of the exact frequency results for the imiform and 
"linear" hinged beams is presented in figure 4 for the case of zero off- 
set. The difference between the resiilts for the two beams is very marked, 
particularly for the first mode. One of the most important things to be 
noted in this comparison is the large difference in slope between the two 
curves for the first mode. The average slope of each of these lines is 
directly proportional to the Southwell coefficient for the first mode 
(see eq. (lO)), and the large difference in slope indicates that a single 
value of this coefficient could not adequately predict the first-mode- 
frequency variations for both beams . This result contradicts the often 
made assumption that the Southwell coefficient is largely independent of 
beam mass and stiffness distribution. 

For the higher modes the slope of each of the lines (fig. 4) is also 
proportional to the Southwell coefficient, but unfortunately each beam 
has a different constant of proportionality. Thus, it cannot be observed 
directly from this figure that the Southwell coefficient for the higher 
modes also varies appreciably with beam characteristics; this fact, how- 
ever, is evident from the charts for frequency determination to be pre- 
sented subsequently. 

Cantilever beams without tip mass .- The frequency of rotating canti- 
lever beams as well as of hinged beams is of interest in the an alysis of 
a teetering rotor because both symmetrical (cantilever) modes and anti- 
symmetrical (hinged) modes may be excited. Consequently, in the following 
paragraphs the Rayleigh approach employed in the present report is eval- 
uated for cantilever beams . 

Frequency resiolts for uniform cantilever beams are presented in 
fig\ire 5. The Rayleigh results are in good agreement with the more exact 
results for the second and third modes. For the first mode, however, the 
maximum error is somewhat larger, about 5 percent in the frequency. 
Nevertheless, the effect of offset on the frequency variation is pre- 
dicted fairly accurately for all three modes. 

For comparison, the results of approximating the first cantilever 
mode by the pendulum mode of a hinged beam are also given in figure 5. 
Frequency results based on this shape are seen to be always less than the 
exact values. As the rotational -speed parameter increases, these res\ilts 
become more and more accurate; for the lower rotational speeds, however, 
the use of the nonrotating -beam first mode shape yields the most accurate 
results. The effects of root offset on frequency are predicted by the 
use of either the pendulum mode or the first cantilever bending mode. 
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The variation of tending frequency with rotational speed is shown 
in figure 6 for a cantilever beam with linear mass and stiffness distri- 
bution and with offsets of 0 and 10 percent. As is the case for the uni- 
form cantilever, the Rayleigh frequency results, based on the nonrotating- 
beam cantilever mode, are very accurate for the second and third modes, 
but are not so accurate for the first mode at the higher values of the 
rotational-speed parameter; however, the effect of the offset is again 
predicted fairly accurately. 

The Rayleigh results based on a penduliom mode, which are also shown 
in figure 6, are again seen to be always less than the exact values and 
to increase in accuracy as the rotational -speed parameter increases. At 
the lower rotational speeds, however, these results are again appreciably 
less accurate than those based on the first cantilever bending mode shape. 
As was the case for the uniform beam, both the pendulum mode and the first 
cantilever mode predict the effects of the offset equally well. 

A comparison of the frequency results for the uniform and "linear” 
cantilever beams with zero offset is given in figure 7 • From the figure 
it is evident that there is only a small difference in the slope of the 
exact first-mode frequency curves and thus in the Southwell coefficient 
for the two beams. This small difference, however, is predicted, although 
not too accurately, by the Rayleigh approach based on the nonrotating -beam 
mode shape; whereas, if a pendulum-mode approximation had been used, no 
difference could have been predicted. 

For the higher modes, the effects of mass and stiffness distribution 
on frequency are more pronounced and lead again to the conclusion that, 
in general, a single value of the Southwell coefficient cannot accurately 
predict the frequency variations for beams with appreciably different mass 
and stiffness distributions. 

The error in the first -mode -frequency results obtained by the Ray- 
leigh approach (fig- 7) is almost the same for both beams. Thus, this 
error apparently is independent of beam. mass and stiffness distribution; 
this observation suggests the possibility of applying a correction, 
based on the known errors for these particular beams, to the Rayleigh 
results obtained for cantilever beams with other mass and stiffness 
distributions . 

Cantilever beams with tip mass .- For beams with a mass at the tip, 

the results for the cantilever case suggest certain simplifications which 
may be carried over to the hinged beams; thus the cantilever results are 
discussed first. 

The variation of bending frequency with rotational speed for a uni- 
form cantilever beam with a concentrated mass at the tip and zero offset 
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is given in figure 8. Results are presented for two cases: tip mass 

equal to the beam mass and tip mass equal to one-half the beam mass. 
Figure 8 shows that the Rayleigh results are of the same order of accu- 
racy as for the beam without tip mass - very accurate for the second and 
third modes but relatively less accurate for the first mode. 


It is of interest to note that for each mode the variation of the 
frequency parameter with the rotational-speed parameter is almost identi- 
cal for the two values of tip mass considered. In fact, if these results 
are compared with those for the beam with zero tip mass in figure 5, the 
variation for all three cases is seen to be practically identical. 


The foregoing observations create the impression that the zero- 
offset Southwell coefficient for each mode is independent of the value 
of the tip mass. This assumption is true for the first mode but is mis- 
leading for the higher modes as is evident from equation (lO) where it 


can be seen that a constant of proportionality 



which varies 


with tip mass, is involved. Nonetheless, inasmuch as this constant of 
proportionality is defined by a ratio of nonrotating-beam frequencies, a 
new rotating -beam frequency coefficient, or modified Southwell coeffi- 
cient Ko ' can be defined (see eq. (ll)) which is essentially indepen- 
^n 

dent of tip mass and, as will be shown subsequently, of beam mass distri- 
bution as well. 


Hinged beams with tip mass .- The variation of bending frequency with 
rotational speed for a uniform hinged beam with a concentrated mass at 
its tip and zero offset is given in figirre 9. Results are given for two 
cases: tip mass equal to beam mass and tip mass equal to one-half the 
beam mass . The Rayleigh results are very accurate for all three modes 
over the entire range of variables investigated, and it may be inferred, 
particularly for the first mode, that the Rayleigh procedure will yield 
reasonably accurate results for appreciably larger values of the 
rotational-speed parameter and tip mass. 

From figure 9 the frequency variation can readily be seen to be con- 
siderably different for the two values of tip mass, unlike the cantilever 
results of figiore 8, for which the frequency variation is essentially 
independent of tip mass. In an attempt to explain this difference between 
the two cases, the resiilts of figure 9 were replotted in figure 10 as a 
function of the rotational-speed parameter used for the cantilever cases, 
that is, the squared ratio of rotational speed to the bending frequency 
of the beam in the first cantilever mode. From figure 10 the frequency 
variation with this rotational-speed parameter may be seen to be essen- 
tially independent of tip mass, as was noted for the cantilever. Conse- 
quently, a new constant which is insensitive to the mass distribution of 
the beam is suggested. For hinged beams this constant is also defined by 
equation (ll) . The invariance of this constant with beam mass distribu- 
tion is discussed subsequently in this report. 
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Estimation of fundamental frequency of beam with tip mass .- A method 
which permits the accirrate approximation of the first bending mode shape 
of a hinged beam with a tip mass from a knowledge of the first mode shape 
of the beam without a tip mass is presented in appendix B. Once such a 
shape has been determined, the computation of the nonrotating-beam first- 
mode frequency and the associated Southwell coefficient is a relatively 
simple matter. In order to illustrate the accuracy of this procedure, 
nonrotating-beam bending frequencies and Southwell coefficients were com- 
puted for the uniform beam with two values of tip mass and were compared 
with the values obtained by using the exact nonrotating-beam bending mode 
shapes . 

For the case of a uniform beam with tip mass equal to one-half the 
beam mass, the nonrotating-beam frequency squared obtained by using the 
approximate shape was found to be too high by about 2 percent, and the 
associated Southwell coefficient was found to be too low by about 2 per- 
cent. If these errors had both been in the same direction, the first 
bending frequency of the rotating beam would have differed by only about 
1 percent or less from the Rayleigh resixlt based on the exact nonrotating- 
beam mode shape. However, because the two errors tend to cancel, the 
difference would be much less. 

The results for the case of a tip mass equal to beam mass showed 
very similar characteristics, although the error in nonrotating-beam 
frequency was slightly higher. 

Although the method of appendix B has been evaluated only for the 
case of uniform beams, it is believed that the method will be equally 
accurate for beams with other mass or stiffness distributions. 


CHARTS FOR BENDING-FREQUENCY DETERMINATION 


In the preceding section, the Rayleigh approach was evaluated and 
the conclusion was reached that Southwell coefficients obtained by using 
nonrotating-beam mode shapes lead to reasonably accurate bending fre- 
quencies of rotating beams, at least for the range of the rotational- 
speed parameter encountered in helicopter blades. The evaluation also 
showed that the Southwell coefficients can vary appreciably with beam 
characteristics . This section describes a group of charts based on the 
Rayleigh approach which permit the rapid estimation of bending frequencies 
of rotating aind nonrotating beams . 
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Rotating and Nonrotating Beams Without Tip Mass 


In order to provide a means for the rapid, yet reasonably accurate, 
estimation of rotor -blade bending frequencies , nonrotating -beam frequency 
coefficients, zero-offset Southwell coefficients, and offset-correction 
factors for the Southwell coefficients have been computed for a series of 
beams with linear mass and stiffness distributions and have been compiled 
in chart form. The range of mass and stiffness distributions was selected 
to encompass variations found in currently manufactured blades with some 
latitude for new design. All the constants are based on the mode shapes of 
the nonrotating beam, which were obtained by standard numerical-iteration 
procedures. (See section entitled "Results for Bending Modes" for more 
details regarding these procedures . ) 

The form of the Rayleigh energy equation which is used in conjunction 
with the charts to obtain bending frequencies is eqmtlon (8) with Kq^, 

K^ , and e in nondime ns ional form: 

-Ln 


where K.. = and K. = K. L. The charts for frequency determination 

On On ^n -^n 

are presented in figures 11 to l6. In each chart, the abscissa is the 
ratio of the beam mass per unit length at the tip of the beam to the mass 
per unit length at the root; 1.0 represents a constant-mass beam and 0 a 
beam in which the mass varies linearly to zero at the tip. Curves are 
presented for three different stiffness variations: the solid curves for 

beams with constant stiffness, the long-dash ciirves for beams where the 
stiffness drops linearly to half the root value at the tip, and the long- 
dash, short-dash curves for beams which have zero stiffness at the tip. 

Each of these curves is faired through only three points, one at 
each end and one at the middle; for the Southwell coefficients and offset- 
correction factors, this procedure should involve little error because, in 
most cases, the variation is nearly linear, but for the frequency coeffi- 
cients the fairing may appear to be questionable. However, the fairing of 
these curves was not entirely arbitrary. The fundamental bending frequency 
of cantilever beams with linear mass and stiffness distributions is given 
in reference l6 for cases in which the mass and stiffness variations are 
proportional^ that is^ where 



( 12 ) 
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Two of the cases considered in this reference, namely, the ones where 
both ratios equal 0 and 1, are identical to cases treated in this report 
and the results for these are in good agreement. The other cases treated 
in reference l6, namely, those for which this ratio is 0.2, 0.4, and 0.6, 
were used in fairing the curves of a^^ for the cantilever case. The 

other curves for the frequency coefficient for cantilever and hinged beams 
were then faired by using this first set of ciurves as a guide. 

Charts which permit the rapid estimation of nonrotating-beam fre- 
quency coefficients, zero-offset Southwell coefficients, and offset- 
correction factors for the Southwell coefficients are presented in fig- 
ures 11 to 13 for beams hinged at the root and in figures l4 to l6 for 
beams fixed at the root. 

Since the zero-offset Southwell coefficient for the pendulum mode is 
always unity regardless of the mass eind stiffness distribution of the beam, 
it is not included in figure 12. However, the offset-correction factor 
for this mode is not independent of mass distribution but is independent 
of stiffness distribution, as indicated in figure 15 • The pendulum-mode 
results in figure 13 are also given in reference 4. 

As was mentioned in the section of this report entitled "Eval\iation 
of Rayleigh Approach, " the zero-offset Southwell coefficients for the 
first cantilever mode (given in fig. 15) will yield accurate rotational 
frequencies only at relatively low values of the rotational-speed parameter 
and must be corrected in accordance with the results of figure 5 or 6 at 
higher values of this parameter. A fixed-percentage correction cannot be 
given because the error is a function of the rotational-speed parameter. 

The effect of root fixity on the Southwell coefficients can be noted 
by comparing the curves of figure 12 with those of figure 15- The first- 
mode results for the cantilever beams shoiild be compared with the pendulum- 
mode Southwell coefficient for the hinged beam which is always unity for 
the case of zero offset. Likewise, the second-mode curves of figure 15 
shoiild be compared with the first-mode curves of figiire 12, and so forth 
for the higher modes . From this comparison it is seen that the effects 
of root fixity on the Southwell coefficients are fairly small and can 
probably be neglected for rough approximations in all cases, except for 
the first cantilever mode. With this assun^ition, the results of figure 12 
for the third bending mode can be used as reasonable approximations for the 
fourth cantilever mode . 

The variation of the Southwell coefficient may be seen from figures 12, 
13, 15, and l6 to be relatively Insensitive to beam stiffness distribution, 
particularly for cantilever beams but also for the hinged beams. This 
observation, coupled with the facts that frequency is proportional to the 
square root of the Southwell coefficient and that the influence of the 
Southwell coefficient decreases for higher modes (for constant rotational 
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speed), leads to the conclusion that fairly good approximations to the 
Southwell coefficients for beams with other than linear stiffness distri- 
butions may also be obtainable from this set of charts. The examples 
presented in the following section appear to bear out this conclusion. 

Application of charts to several actual helicopter blades .- To illus- 
trate the use and the type of accuracy which can be expected from the 
frequency charts of figures 11 to l6 and to demonstrate that the charts 
work well even when the mass and stiffness distributions of the beams are 
not exactly linear, bending frequencies have been estimated for the first 
three modes of four existing helicopter blades, all of which are hinged. 

The following procedirre, which may be made clearer by reference to the 
sketches in table I, was used in the estimation: 

(a) Straight lines were faired throiigh the mass and stiffness dis- 
tributions for the blade; large values near the root were ignored. 

(b) From these fairings, the effective root values mQ and EI^ 
and the necessary tip-root ratios were obtained. 

(c) By using these ratios, values of a^^, Kq, and K-]_ were 
obtained from the appropriate charts (figs. 11 to l6). 

(d) Substitution of these constants and e into the Rayleigh equa- 
tion (eq. (12)) yielded the bending frequencies at zero and the rated rotor 
speed. 

The mass and stiffness distributions for the blades considered are 
shown on the left side of table I. The actual distribution is given by 
the solid lines, and the linear approximation, selected to represent this 
variation, is given by the dashed lines. These linear approximations 
used in estimating the frequencies were the initial ones selected, and 
no attempt was made to improve them in order to obtain the best agreement 
for all modes. The frequencies shown as "exact" in table I are values 
furnished by the manufacturer. 

A comparison of the exact and estimated results given in table I for 
these blades indicates that the estimated results are very accurate when 
the crudeness of the linear approximations used is considered. 

Although no comparisons have been made for cantilever blades because 
sufficient information regarding such blades was not available, even more 
accurate results should be obtainable for this end condition since large 
values of root stiffness can be taken into account more accurately by con- 
sidering the blade to be cantilevered at the outboard edge of the stiff 
region and then using the offset-correction factor for the Southwell 
coefficients . 
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Beams With a Mass at the Tip 

Uniform cantilever beam .- Expressions defining the bending frequencies 
and mode shapes of nonrotating uniform cantilever beams with a tip mass 
equal to a fraction r of the beam mass are given in reference I7. These 
expressions^ in somewhat simpler form, are the following: the defining 

relation for the frequencies is 

1 + cos 0 cosh 0 - r0(sin 0 cosh 0 - cos 0 sinh 0) = 0 (13) 


where 






and the mode shapes are 

sinh + sin 0^ 

y (x) = sinh x - sin x + (cos x - cosh x) (l^) 

cosh + cos 9j^ 

In addition to the defining relation for the frequency, reference I7 
also gives values of 0^ for the first three modes of cantilevers and 

for several values of r. Some of these results, which are pertinent to 
helicopter blades, are plotted in figure I7. Values of rather than 

0^ are plotted, because is directly proportional to frequency and 

corresponds to the nonrotating frequency coefficients a^ presented 
previously. 

For larger values of r fairly accurate values of 0^2 can be 
obtained from the following approximate expression: 



(15) 


where is a constant for each mode which can be determined from the 

frequency resiilts for the largest value of r - in this case, 2. Equa- 
tion (15) can also be used for nonuniform beams and for hinged as well as 
cantilever beams. 
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The variation of the Southwell zero-offset coefficient with tip mass 
is given in figure I8 for the first three modes of a uniform cantilever 
beam. These resxilts were computed by using the mode shape of the nonro- 
tating beam presented in equation (l^); the integrations were performed 
analytically. Although only three points were used to establish each 
curve of figiire I8, the fairing shoiild be quite accurate since the vari- 
ations shown are almost linear. 

The Southwell coefficients of figiu-e 18, in conjunction with the 
nonrotating-beam frequency coefficients of figure I7, should permit very 
accurate estimates for the bending frequencies of rotating uniform beams 
with a tip mass except possibly for the first mode, for which a correction 
may be made in accordance with results given in figure 8 for large values 
of the rotational -speed parameter. 

The effect of root offset has not been studied for this case, but 
offset -correction factors can be obtained from the mode shapes defined by 
equation (l4). 

Uniform hinged beam.- By using the method of reference I7, expressions 
defining the bending frequencies and mode shapes of nonrotating uniform 
hinged beams with a mass at the tip have been determined. The defining 
relation for the frequency is 

2r0 + coth 0 - cot 0=0 (16) 


and the mode shapes are given by 


y„(x) 


sinh X + 


sinh 0^ 
sin 0^ 


sin X 


(17) 


Values of 0^^ have been determined for several values of r; these 
resiilts are given in figure 19 as frequency coefficients 0^^^, together 
with the frequency coefficients for the case of zero tip mass. 

By using the nonrotating-beam mode shape, given by equation (l7)j 
values for the zero-offset Southwell coefficient have been determined 
for hinged beams with a tip mass and are given in figiare 20. For the 
pendulum mode, Kq is always imity and therefore is not shown. The 

results in fig\ares 19 and 20 together permit the rapid estimation of the 
bending frequencies of rotating uniform hinged beams with a mass at the 
tip . 
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Pendulum -mode results for hinged beams vith linear mass distribu- 
tions . - The zero-offset Southwell coefficient for the pendulum mode of 

a hinged beam is equal to unity, regardless of the mass or stiffness 
distribution of the beam. For the case of hinge offset, however, the 
Southwell coefficient is independent of stiffness distribution but varies 
considerably with beam mass distribution and with the tip mass. A chairt 
(see fig. 21) has been prepared which permits the rapid estimation of the 
offset-correction factor to the Southwell coefficient for hinged beams 
with an approximately linear mass distribution plus a mass at the tip. 

First bending mode frequency of nonuniform hinged beam . - A simple 

method is indicated in appendix B for obtaining an approximate first 
mode shape for any beam with a tip mass from a knowledge of the beam 
mode shape without a tip mass . Once such a shape is determined, the 
fundamental bending frequencies of the rotating eind nonrotating beams 
can be determined very easily by application of the Rayleigh frequency 
equation (eq. (l)). 


Rotating Beams With Nonlinear Mass Distribution and 
Approximately Linear Stiffness Distribution 

In the section of this report concerned with the evaliaation of the 
Rayleigh approach, a modified form of the zero-offset Southwell coeffi- 
cient shown to be insensitive to variations in beam tip mass. 

This coefficient is defined for both cantilever and hinged beams by equa- 
tion (11). 

In order to determine whether this new coefficient is also insensi- 
tive to other variations in beam mass distribution, all values of Kn 

'-’n 

presented in the charts for rapid frequency estimation were converted to 

Kp, ' . For each stiffness distribution IC-i ' was found to be almost con- 
^n ^n 

stant for each mode , the differences being of the sajne order of magnitude 
as the errors inherent in the Rayleigh approach used herein. 

To facilitate the estimation of bending frequencies for rotating 
beams with large tip masses or possibly other nonlinear mass distribu- 
tions, values of ^" 0 ^ 8,11 the beams treated in the present report 

are plotted in figure 22(a) for cantilever beams and in figure 22(b) for 
hinged beams . Curves have been faired through the points to give average 
values for 8,nd thus for Kq for beams with approximately linear 

stiffness distributions and with any mass distribution. In analyzing 
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these results, the facts that frequency is a function of the square root 
of K„ and that the influence of on frequency decreases with 

increase in mode number should be kept in mind. 

From equation (ll) it is apparent that the first bending frequency 
of the nonrotating beam cantilevered at the root and the nth bending fre- 
quency of the nonrotating beam with its actual end fixity are required 

to determine Kn (and thus the bending frequency of the rotating beam) 
^n 

from a knowledge of Kn ' . In spite of this complication, however, the 

n 

charts presented should be useful in design studies involving rotating 
beams with nonlinear mass distributions but with approximately linear 
stiffness distributions. It should be emphasized at this point that the 
constancy of Kq ' has been demonstrated for only a limited variety of 
^n 

mass distributions, and thus application to blades having mass distribu- 
tions radically different from those considered in this report should be 
made with caution. 


Rotating Beams With Mass and Stiffness Distributions 
Not Representable by Foregoing Approximations 

The charts presented in this report permit the rapid estimation of 
bending frequencies for rotating beeims with a mass and stiffness distri- 
bution each of which cein be reasonably approximated by a straight line and 
for uniform beeims with a tip mass; also the charts facilitate the estima- 
tion of bending frequencies for rotating beams with fairly arbitrary mass 
distributions and approximately linear stiffness distributions. For 
other cases, for example, beams in which the stiffness varies irregularly 
n.1 1 along the blade, the basic Rayleigh energy method utilizing the modes 
of the nonrotating beam may be used. Although this method has been eval- 
uated in this report only for linear distributions of mass and stiffness 
and concentrated tip mass, there is no reason to believe that it will not 
work equally well for other distributions. All that is required in this 
approach is the frequency and mode of the nonrotating beam, which can be 
determined by methods such as are described in references 2 and 15 . (A 
method which gives directly the required first derivative of the mode as 
well as the mode shape itself is preferable . ) With such res\ilts the 
integrals of equation (l) can be evaluated readily by accurate n\jmerical 
methods such as those of reference 15 , and values can be obtained for the 
Southwell coefficient from which the bending frequencies at any rota- 
tional speed can be determined with little effort. 
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Mode -Expans ion Method 


A more acciorate mode -expansion method for determining the bending 
frequencies and modes of a rotating or nonrotating beam has been devel- 
oped in appendix B and has been used as a yardstick in the evaluation of 
the Rayleigh approach. In this approach the lowest three bending modes 
and frequencies are obtained by the solution of a fifth-order determlnantal 
equation for cantilever beams and a sixth-order equation for hinged beams. 
In order to facilitate the further application of this method to the 
accurate determination of the modes and frequencies of rotating and non- 
rotating beams, certain integrals which have been evalijated sire presented 
in table II. These results permit the setting up of frequency determi- 
nants for beams with any combination of linear mass and stiffness distri- 
bution, concentrated tip mass, offset, and rotational speed (including 
many combinations not treated herein) . With the evaluation of additional 
integrals (some of which are given in ref. l8), these results can be used 
to determine the bending frequencies and modes for rotating and nonrotating 
beams with concentrated mass at other locations or with higher order mass 
and stiffness distributions. If practice dictates the necessity of addi- 
tional charts for other combinations of linear mass and stiffness distri- 
bution and tip mass or for parabolic beam mass and stiffness distributions, 
it might be advantageous to use this method to set up such charts if high- 
speed computing machines suitable for solving the determlnantal eqiiations 
are available. 


The frequency charts and procedures for frequency determination of 
this report have all been directed toward the determination of frequencies 
for uncoupled bending vibrations perpendicular to the plane of rotation. 

In cases where the principal axis of the blade cross sections (axis about 
which the stiffness is a mlnimimi) is not parallel to the plane of rotation, 
natural bending vibrations having the lowest frequency will take place 
perpendicular to the chord. An extreme case of such vibrations would 
occur if the blade chord were perpendicular to the plane of rotation, in 
which case, blade vibrations woiild take place in the plane of rotation. 

Frequencies of vibration, when the blade chord is Inclined at any 
angle with the plane of rotation, can be determined from the fre- 
quencies of vibration perpendiciilar to the plane of rotation by means of 
a simple formula proposed in reference 19: namely. 


Vibration in Planes Other Than Those Perpendicular 


to Plane of Rotation 
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where is the frequency of bending vibrations perpendicular to the 

L 

plane of rotation and is the frequency for bending vibrations in a 

plane making an angle with the axis of rotation. 


At large angles of attack, the indicated correction may be signifi- 
cant for the lower modes. However, inasmuch as ^ is usmlly 5 "to 


10 times as large as for the lowest bending mode of helicopter blades 

and even larger for the higher modes, in most cases the angle of attack 
of the blade will have little effect on bending frequency and may be dis- 
regarded. This fact is significant since it indicates that blade fre- 
quency will not change appreciably during each revolution because of 
cyclic-pitch changes and thus may be assumed to be constant. 


RESULTS FOR BENDING MODES 


In the process of obtaining the frequency results presented in the 
preceding sections of this paper, a large number of mode shapes of both 
rotating and nonrotating beams with various mass and stiffness distribu- 
tions were determined. These results are presented in tabular form in 
order to make them more useful in analytical studies and are compared in 
this section with each other in order to show the effect of the various 
parameters on mode shape . 


Nonrotating Beams 

The first three mode shapes for nine nonrotating cantilever and nine 
nonrotating hinged beams with different combinations of linear mass and 
stiffness distributions are given in tables III and IV, together with their 
first and second derivatives. These results were obtained by standard 
numerical-iteration procedirres. For the cantilever beams (table III), the 
procedure of reference I 5 was used with 10 stations; step-integration 
procedures were used for the first mode, and equivalent -load methods were 
used for the second and third modes. For the hinged beams (table IV), a 
matrix-iteration procedure using weighted integration matrices similar to 
those given in reference 21 was employed with I 5 stations. More stations 
were needed for the hinged beams than for the cantilever beams because the 
third hinged mode has one more loop or node than the third cantilever mode. 

In order to illustrate the accuracy of the nonrotating mode shapes 
computed by this method, the exact results given for the uniform beam in 
reference 20 are also included in tables III and IV. A comparison of the 
results indicates that the error of the present results is less than 
1 percent. Nonrotating mode shapes are shown for the hinged beams in 
figure 25 and for the cantilever beams in figure 24. 
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Rotating Beams 

The mode and frequency results for rotating beams were obtained in 
the present paper by the method of appendix A. This method yields mode 
coefficients which, when multiplied by the mode shapes of nonrotating 
uniform beams normalized to positive tip values and summed, give the mode 
shapes of the rotating beam. These coefficients can also be used in con- 
junction with the spanwise derivatives of the uniform-beam mode shapes to 
obtain similar derivatives for the rotating beams . The required uniform- 
beam modes and derivatives are given in reference 20, but they are not all 
normalized to positive tip deflections and thus certain sign modifications 
are necessary. These modes and the first two derivatives are also given 
in tables III and IV with the proper signs and tip deflections . 

All the mode coefficients for rotating beams obtained in the present 
investigation are given in tables V and VI. These coefficients have been 
normalized in such a manner that the modes obtained by using them will 
have the same tip deflection as the uniform-beam modes used in the compu- 
tation. Table V contains the resiUts for the hinged beams, whereas 
table VI contains those for the cantilever beams. 


Comparison of Rotating and Nonrotating Beams 


Hinged beams .- The mode shapes of a uniform hinged beam for zero 
rotational speed and a rotational speed Q equal to the first bending 
frequency are shown in figure 25. A comparison of these shapes 


indicates that although some differences between the modes exist, they 
are relatively small, particularly for the higher modes. 


A similar comparison is given in figure 26 for hinged beams with 
linear mass and stiffness distributions, both zero at the tip. For this 
case the difference in mode shapes is very small for all three modes; 
this undoubtedly accounts for the fact that the Rayleigh approach was 
found to be very accurate for this case. (See fig. 3-) 


By comparing the results of figiires 25 and 26, a large disagreement 
may be noted between the mode shapes of the two beams; this disparity 
apparently accounts for the substantial differences in the Southwell 
coefficients for the two beams . 


The calculated mode shapes have not been plotted in a form which 
shows the effect of offset on the mode shapes of rotating beams; but by 
comparing the mode coefficients for 0- and 10-percent offsets in 
table V, the effect may be seen to be small. 
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Cantilever beams.- The modes of rotating and nonrotating uniform 
cantilever beams are shown in figure 27- From the figure the mode shapes, 
particularly those for the first and second modes, may be seen to change 
appreciably with rotational speed. 

A similar comparison can be made for cantilever beams with linear 
mass and stiffness distributions on the basis of the results shown in 
figure 28. The mode shapes vary in about the same marmer with rotational 
speed for this type of beam as for the uniform beam. 

If the mode coefficients for 0- and 10-percent offsets in table VI 
are compared, the effect of offset on mode shape is again seen to be very 
small for both beams. 

Beams with a mass at tip.- Bending mode shapes for a rotating and a 
nonrotating uniform hinged beam with a mass at the tip eqml to the beam 
mass are shown in figure 29. The differences in mode shape are very 
small for all three modes. This similarity apparently accoxmts for the 
excellent accuracy of the Rayleigh approach for this configuration. 

Similar results for a uniform cantilever beam with a mass at the tip 
equal to the beam mass are presented in figure 30. For this case, results 
are given for three values of the rotational-speed parameter, namely, 

^ = 0, 10.^5, and lk.^6, and also for the nonrotating viniform beam 

without tip mass. From this figiire the rotating -beam mode shapes may be 
seen to be only slightly different from each other but considerably dif- 
ferent from the nonrotating shape, particularly for the first and second 
modes, and vastly different from the mode shape of the beam without a 
tip mass. 

Mode coefficients for rotating uniform hinged and cantilever beams 
with a mass at the tip are listed in tables V and VI. Mode shapes for 
nonrotating uniform beams with a mass at the tip have not been tabulated 
but can be calculated by means of equations (l4) and ( 17 ) for any value 
of tip mass . 


CONCLUDING REMARKS 


A Rayleigh energy approach, which utilizes the mode shape of the 
nonrotating beam as an approximation for the mode shape of the rotating 
beam in the determination of the bending frequencies of the rotating 
beam, has been evaluated. The evaluation led to the conclusion that 
this approach yields reasonably accurate bending frequencies for rotating 
hinged and cantilever beams with arbitrary stiffness and mass distribu- 
tions, including concentrated masses, at least within the limits of the 
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rotational speeds currently encountered by helicopter blades. The evalu- 
ation also showed that the Southwell coefficients vary appreciably with 
beam mass distribution and, to a less extent, with beam stiffness distri- 
bution. A modified form of the zero-offset Southwell coefficient, which 
involves the nonrotating -beam frequencies, was found to be insensitive to 
changes in beam mass distribution. 

By using the Rayleigh approach as a basis, several groups of chants 
and associated procedures have been presented, which permit the rapid 
estimation of the first three bending frequencies for a variety of 
rotating and nonrotating hinged and cantilever beams . Since the charts 
are not applicable to all beams, practice may dictate the need for addi- 
tional charts which may be set up by vising the methods described. The 
charts and associated procedures presented in this report are summarized 
below, the most easily applied being listed first; 

(a) Charts are presented which permit the rapid estimation of bending 
frequencies of rotating and nonrotating beams with mass and stiffness dis- 
tributions, each of which csui be approximated by a linear relation. In 
example applications, this procedure has been shown to give good results 
for the bending frequencies of several actual helicopter blades with mass 
and stiffness distributions appreciably different from linear. 

(b) Charts are presented for rapidly estimating the effects of tip 
mass on the rotating and nonrotating bending frequencies of uniform beams. 

(c) A chart is presented which permits the rapid estimation of the 
effects of offset on the pendulum frequency of hinged beams with any stiff- 
ness distribution, an approximately linear mass distribution, and a con- 
centrated tip mass. 

(d) A simplified procedure is presented for estimating the first 
bending mode and frequency of a rotating or nonrotating hinged beam with 
a tip mass from a knowledge of the first mode shape of the nonrotating 
beam without a tip mass. 

(e) Charts for a modified Southwell coefficient, which appears to be 
insensitive to changes in beam mass distribution, are presented; these 
charts permit the rapid estimation of the first three bending frequencies 
of rotating beams with approximately linear stiffness distributions from 
a knowledge of the bending frequencies of the nonrotating beam. 

(f) Bending frequencies for beams with unusual mass and stiffness 
distributions which cannot be estimated by using the charts can be deter- 
mined directly from the Rayleigh energy equation by first calculating the 
bending frequencies and associated mode shapes of the nonrotating beams. 
This approach can be expected to yield results which are in error by less 
(usually much less) than 5 percent, except for the first cantilever fre- 
quency which may be in error by as much as 5 percent but which can easily 


26 


NACA TN 5^59 


be corrected to give a much more accurate result. The method has the 
adAraintage over other simplified approaches of improved accuracy and 
wider applicability and over more exact approaches of simplicity and 
flexibility. 

A more accurate mode -expans ion method for determining the bending 
frequencies and modes of a rotating or a nonrotating beam has been devel- 
oped and has been used to evaluate the Rayleigh approach. In order to 
facilitate the f\orther application of this method to the accurate deter- 
mination of modes and frequencies of rotating and nonrotating beams with 
combinations of linear mass and stiffness distribution and concentrated 
tip mass different from those considered herein, certain integrals which 
have been evaluated are presented in tabular form. 

In conjunction with obtaining the frequency resixLts which comprise 
the greater part of this report, bending mode shapes were determined for 
a wide variety of hinged and cantilever beams. These results show the 
effect of rotational speed, mass and stiffness distributions, offset, 
root fixity, and other parameters on bending mode shape; they have been 
tabulated in normalized form together with their first and second deriva- 
tives or as mode coefficients which, in conjunction with tabulated modes 
and derivatives of uniform beams, permit the rapid determination of the 
mode shape and higher derivatives as well. The tabulated results should 
prove useful in other analyses, for example, in the simplified approach 
presented in an appendix. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February 2k, 1955 - 
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APPENDIX A 

SOLUTION OF DIFFERENTIAL EQUATION FOR ROTATING BEAM 
BY EXPANSION IN TERMS OF NORMAL MODES OF 
UNIFORyl NONROTATING BEAM 
Solution by Galerkin Method 


The eqmtion of motion which defines the bending vibrations perpen- 
dicular to the plane of rotation of a rotating beam with a concentrated 
mass at its tip can be written as 





(Al) 


where 


and 


8(x-L) =0 (x ^ L) 


6(x-L) = - (x = L) 
L 


T = 



e)m dll + Mi;(L + e) 


or, in nondimens ional form. 



t»n^yn + rbn%(x-l)yn(l) 



(A2) 
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where 


6(x-l) =0 (x / 1) 

5(x-l) =1 (x = l) 


and 


Ti = 



(t] + e)m(Tj)dfi + r(l + e) 



EIo 


Each normal mode of the rotating beam can be expanded in te^ of 
the modes of a uniform nonrotating beam with the same end restraints 
as follows : 

00 

yn = ^ (^5) 

q=0 ^ 

where the quantities axe the normalized bending mode shapes of a 

stationary uniform beam, and the coefficients are ^determined. 

Substituting this expansion into equation (A2) gives 


dx^ 


El ^ 

dx2\q=0 


qrq 


q=0 q=0 ^ 




= 0 


(a4) 
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One way of determining the coefficients A^^ from this equation 

is the Galerkin procedure which consists in mxiltiplying the equation 
hy and integrating over the length of the beam. Thus, 



(A5) 

Integrating the first term in equation (A5) by parts twice and the 
last tern by parts once and making use of the known boundary conditions 
gives for either a cantilever or a hinged beam: 

/ 9*p"S Zl A. 62 f S9!p X 62 - 

^ 0 q=0 ^ ^ 0 q=0 ^ 


rbn%(l)Z. An 9<q(l) + (^\ f A ?(q' dSE = 0 (a6) 

^ q=0 ^ \^1 y '^0 q=0 ^ 

where the primes designate differentiations with respect to x. Inter- 
changing the order of integration and summation yields; 





rbn%(l)2l Anq0q(l) + 

q=0 





0 




dx = 0 


(A7) 
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Equation (A7) can be rewritten as 




q=0 


Iqp " V Mqp 


+ Rqp) + 


V%Ri, 




= 0 


(A8) 


in tenns of a new set of constants: namely, 

Iqp = r S^p'Vq" 
^0 



- i’9^p(l)9^q(l) 


(A9) 


Sqp - 



Tp^^pVq' ^ 


These coefficients are symmetric; that is, Iqp = Ipq, and so forth. 


For practical purposes^ the expansion must be limited to a finite 
number of nonrotating uniform-beam modes. In this case the summation 
goes from q = 0 to m and equation (a 8) yields m + 1 equations of 
the form 


w 

An^Bqp = 0 

q=0 


where 


®qp = ^qp “ ^n (^qp + ^qp ^ + 


ft 




^l^Sqp 


so that the coefficients Bqp are also symmetric 
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The modes and frequencies of the system represented by this group 
of equations can be obtained for any value of the rotational-speed param- 
eter by equating the following determinant of the multipliers B 

of the mode coefficients Aji^ to zero; 


Boo 

BOI 

Bo2 • • • 

BOm 

®10 

Bll 

®12 • • • 

Blm 

®20 

®21 

^22 • • • 

®2m 


= 0 (All) 


%0 


®ml ®m2 • • • ®mm 


This determinantal equation can be solved by trial and error, with 
any method of evaluating determinants, such as Grout's, to obtain the 
frequency coefficients bj^ and subsequently the associated mode coef- 
ficients An^ for a rotating beam. The resonant frequencies for Ip, 

2p, or np resonant conditions can also be obtained directly from the 
determinant. For small values of than about 0.8, solutions 

can also be obtained by the matrix- Iteration procedure; for larger values 
however, convergence is poor, and \indesired negative values of the fre- 
quency squared (imaginary frequencies) may be encountered before the 
desired positive values are obtained. In the present investigation 
the frequency determinants (eq. (All)) were solved by trial -and-error 
methods with automatic computing machines of the punchcard type. 

For the case of a beam without a tip mass, r = 0, and thus Rgp 
is not needed and Sgp is simplified slightly. If, in addition, the 
beam is miiform, Iqp and Mgp are zero by orthogonality for q ^ p; 

thus for this case the lonknown frequency coefficients bn occur only 
on the principal diagonal. If the determinantal equation is divided 
by then for this case the rotational-speed parameter also 

appears only in the terms on the principal diagonal. 


52 


NACA TN 3^59 


Evaluation of the Integrals Iqp, Mgp, Rgp, and Sqp 

The integrals Igp, Mgp, Rqp, and Sqp may be evaluated mmerically 

by a method such as that given in reference I5, or, if the mass and stiff- 
ness distributions of the beams are defined by analytical expressions, 
they cein sometimes be evaluated in closed form. (See ref. 1, pp. 533“536> 
for instance . ) In some cases integrals already evaluated and tabxilated 
in reference 18 can be employed; these results, converted to the coor- 
dinate system and tip deflection of the present paper, were employed where- 
ever possible in the present study. In this report a3JL integrals for 
the uniform rotating beams with and without a tip mass were evaluated by 
exact methods. Some were also evaluated by numerical methods in order 
to determine how many stations were required to obtain good accuracy. By 
this procedure about 25 stations were found to be required for some of 
the integrals involving the fourth and fifth modes. 

For the nonuniform rotating beams, Iqp, Mqp, and Rqp were evalu- 
ated both exactly and nmerically, but Sqp was evaluated only numeri- 
cally because of the effort involved in evaluating this integral exactly. 
All the integrations performed in this report are based on mode shapes 
normalized to unity at the tip. Where numerical integrations were made, 
the mode shapes and derivatives were obtained from reference 20, but the 
results were modified to correspond to shapes with a unit positive tip 
deflection. 

The remainder of this appendix is devoted to the presentation of 
results (in both nmerical and analytical form) for Iqp> Mqp, I^qp> 

and Sqp which were obtained in connection with the present study but 

which are also applicable to cases not treated in this report. 

Numerical res-ults for beams with linear mass and stiffness distri- 
butions with or without tip mass and offset .- If only linear variations 
in beam mass and stiffness are considered and if they are expressed as 


m = niQ(l - kx) 

El = EIo(l - c5f) 


(A12) 


then the various integrals can be evalmted expeditiously by splitting 
them up as indicated in the following equations: 
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Spq = Spq^ - kSpqj^, + e^Spq^ - kSpq^e) + + ®)Spq-t 


T = T - cT 

-^qp -^qPo qPc 


Mqpj^ - Mqp^ - kMqpj^ 
Rqp = r 


(A13) 


where, in txirn. 



(Al4) 
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All these integrals are obviously symmetric in p and q. Numerical 
values for them are given in table II for values of p and q from 
0 to 5 for hinged beams and 1 to 5 for cantilever beams. As may be seen 
from equations (AI3), these results permit the rapid calculation of the 
terms of a frequency determinant for a rotating beam with any combina- 
tion of the following parameters; (a) linear mass distribution, 

(b) linear stiffness distribution, (c) any offset (including large 
values), and (d) any tip mass. In addition, the resialts can be used 
in conjunction with values of additional integrals to set up similar 
determinants for beams with higher order mass and stiffness distribu- 
tions and beams with concentrated masses at other locations. 

Integrals for \miform beams with tip mass .- In order, to facilitate 
the extension of the results for the ttnlform rotating beams to higher 
modes, the exact expressions for integrals pertinent to such cases are 
included herein. 

The integrals for the cases where p = q can also be used to 
determine values for Southwell coefficients for modes higher than the 
third. The integrals were evaluated by the method of reference 1 or 
taken from reference I8 and transformed into the notation of this report. 
The expressions are given in terms of the parameters 03, 7s> 

values of the first two can be obtained from reference 20 for values 
of s from 1 to 5. For s > 5, Og = 1 for all practical purposes 

and 3s can be obtained from the appropriate frequency equation for the 
nonrotating uniform beam. The sqiiare of 3s the frequency coeffi- 
cient for the nonrotating beam as for the sth bending mode of a uniform 
beam. Values of 7s required for the cantilever beams; for^ 

hinged beams, 73 = 1 for s > 3j the values for s 3 are given in 
the following table; 


s 

7s 

1 

1.02827 

2 

1.00121 

3 

1.00005 
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Tension integrals Sgp for cantilever beams: 

If V 4 <1} 


Sqp - 


2 2 
Pq Pp 

Pq - Pp 


l(.l)'^(=p|3p - OqPq) - 


(- 1 ) 1*1 


apPp - OqPq + 


Kp < 




h 4 
Pq - Pp 


^Pq^Pp^ 

P ^ - P ^ 
Pq Pp 


r(l + e) 



•jopPp 


3^^ + (.l)P^p2pq2 


- °^Pq 




If p = q, 


Sqq 


aqPq 













1 

— + 
2 
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Tension integrals Sqp for hinged heams: 
If p / q / 0, 


4P(i^Pp^ 

r« 

e 

‘*Sq‘‘Pp‘‘ 


-=!■ 

l 

CQ. 

Pq'^ - Pp'' 



(-DP+'l 


PqPp^yp i 

\ 



+ 


r(l + e) 






If p ^ q, hut p = 0, 


Sao = e(-l)^ + r(l + 5) 

Pq# 


If p = 


Sqq - 


^.4 



If p = q = 0, 


Soo =— +— +r(l+e) 
uu j 2 
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Stiffness integrals Iqp for cantilever or hinged beams: 

If q ^ P, 


Iqp = 0 


If P = q 0, 



If p = q = 0, 


IqO = 0 


Mass integrals Mgp for cantilever or hinged beams; 


If V 4 1) 


Mqp - 0 


If P = q ^ 0, 



If P = q = 0, 

Moo * 5 


Tip-mass integrals for cantilever or hinged beams: 
If p = q, or p 7^ q, 


^qp ~ 


r 
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APPENDIX B 


AN APPROXIMATE METHOD OF OBIAINING FIRST BENDING MODE OF 


HINGED BEAM WITH TIP MASS FROM FIRST BENDING MODE OF 


BEAM WITHOUT TIP MASS 


The vibration modes of a rotating hinged beam must satisfy the 
following equation, which expresses the condition of zero moment at the 
root : 


For any given beam the mode shapes of the nonrotating beam can 
readily be shown to satisfy this criterion exactly if e is zero and 
very closely if e is small; therefore, the nonrotating-beam mode shapes 
are good approximations to the rotating -beam mode shapes, regardless of 
the mass distribution of the beam. However, the nonrotating-beam mode 
shape must be that of the beam with the same mass distribution; the pur- 
pose of the present derivation is to go a step further and to obtain an 
approximate first mode shape for a nonrotating beam with tip mass in 
terms of the first mode of the same beam without tip mass. In view of 
the preceding argument, the mode shape obtained in this manner should 
serve as a good approximation to the first mode of the rotating or 
nonrotating beam with the same tip mass and when used in conjunction 
with the Rayliegh approach (eq. (l)) should yield a good approximation 
for the first bending frequency of a rotating or nonrotating hinged beam. 

In deriving such a relation the assumption is made that the second 
derivative or curvature of the beam remains unchanged in the two configu- 
rations. Thus, the mode shape for the beam with tip mass is assumed to 
be of the form 



(Bl) 


or, in dimensionless form 



(B2) 


Yl* « Yl + DqX « + DqX 


(B5) 
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where the first mode of the rotating beam without tip mass is 

assumed to be approximately equal to the nonrotating-beam first mode 
shape Y^. With this mode shape, the criterion of equation (B2) becomes 


2 nl Pi 

/ mx(Yq + Dqx)cLx - / m(x + e)(Y 2 ^ + DQx)dx = 0 (b4) 

^0 0 

If, now, the mass distribution is considered to be made up of the con- 
tinuous distributed mass of the beam m^j^ plus a concentrated tip mass, 

equation (b 4) can be written as 




m<ix(Yi + DQx)dx + 



Yi(l) + 



+ e)(Yi + I\)x)dx - r(l + e) 


Yi(l) + Do 


= 0 


(B5) 


Inasmuch as Yq and x (the pendulum mode shape) are mode shapes of 
the hinged beam with mass distribution mq, they must satisfy the 
orthogonality condition for normal modes for such a beam, namely. 


mqxYq c3x = 0 


and, hence, equation (B5) becomes 




Dd / %x dx + ( — 1 r 


Q. 


'0 


Yq(l) + Do 


-T)q m^x^di? 
'^0 


e f m^Yq dx - eDo r 5^^ dx - r(l + e) 
Jq ^0 


Yq(l) + Do 


= 0 (b6) 
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When this equation is solved for the result is 


= T 


e r m^Yq dx + erYq(l) - 
0 

(?r ■ >1 

rYi(l) 

[(?j - 

- -2 - 

/ m^j^x dx + 

^0 

|(?r - ‘] 

- - - 

r - e / mjj^x dx - er 

'^0 


(B7) 


If the offset e is zero, equation (B7) takes the much simpler form 

-rYi(l) 


Eb = 


(b8) 


0 


m^j^x^dx + r 


or, with nonneilized to unity at the tip. 


Do = 


1 + 


-1 


(B9) 


0 


— P - 

m^x dx 


By comparing the relative values of the terms of eqmtion (B7) and by 
considering the overall influence of terms containing e, small offsets 
can be shown to have a negligible influence on the value of Dq for 

values of the rotational-speed parameter encountered in helicopters. 

Also, for nonrotating beams, e does not enter the problem and, hence, 
can be set equal to zero; thus, as mentioned before, the mode shape, 
based on the result of eqmtion (B9) , obtained in the following paragraphs, 
should serve as a good approximation for both rotating and nonrotating 
beams with and without offset. 

Upon substituting the value of Dq in equation (B9) into equa- 
tion (B3) , the desired first mode shape of the beam with a mass at the 
tip is obtained as 

yi* = s 


(BID) 


_p _ 

m,x (ix 
d 


1 + 
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and the slope J and curvature ( y^* J' of this mode shape are then 

given by 



(Bll) 



(B12) 


(Eq. (b 12), of coiirse, expresses nothing more than the assumed equality 
of the second derivatives.) If the mode shape of a beam with a particialar 
mass and stiffness distribution (but without tip mass) is known, expres- 
sions (bio) to (B12) thus penult the determination of an approximate mode 
shape for the same beam with any concentrated mass at the tip and can be 
used to evaluate the Integrals of the basic Rayleigh equation (eq. (l)) 
by numerical methods j reasonably accurate values can easily be obtained 
in this manner for and for and from which the bending 

frequency at any rotational speed can be determined directly. 


Beams With Linear Mass Distribution Plus Tip Mass 


For the particular case of beams with a linear mass distribution plus 
a tip mass, m^ = 1 - kx and 




x^dx 


1 


3 


k 

k 


Thus 


Ib = 


-1 


1 + 


4 - 5k 


(B13) 


12r 


k2 
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ircsuQ.'fc cstn t )0 usocL in con^^nnciion wiili "tiic firsi nioi.0 sti8.p0 giv0n 
foi* hinged. lD09ins wiili iin0sx mnss sind. stiffness distributions in tnble IV 
to obtain mode and frequency results for such beams . 


Beams With Uniform Mass Distribution Plus Tip Mass 

For the case of beams with a uniform mass distribution plus a tip 
mass, 5(3^ = 1 and thus 


Do = 



(Bl 4 ) 


Uniform Beam With Tip Mass 

For the case of a uniform beam with an arbitrary tip mass. 


% = 1 


and 

El = 1 

Thus, Dq is the same as for the preceding case. In this special case 

the integrals of the Rayleigh equation (eq. (l)), which permit the deter- 
mination of and K and thus of cdr, can be evaluated exactly by 

the methods of reference 1 or I8. The results are 





r 

5r + 1 


(B 15 a) 



El 




(BI5B) 
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- 

I” 

/ 

(^1*) 

^ 0 

V / _ 



12 


+ 






(B15c) 



where = 5*9266, from the results given in reference 20. In the 

preceding integrations On (ref. 20) has been taken equal to unity; 
this assumption introduces a small error of less than 0.1 percent. 


Equations (BI 5 ) are based on Yj rather than normalized to 

unity at the tip. To obtain equivalent formulas for y^^* normalized to 

/ o \ 2 


unity at the tip, these resiilts must be divided by the factor 


3r + 1 


Nonrotating- and rotating-beam frequencies obtained by this method 
for the uniform beam are compared with more accurate results in the sec-^^ 
tion of this paper entitled "Charts for Bending -Frequency Determination." 
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TABLE I 

EXACT AND ESTIMATED FREQUENCIES FOR SEVERAL 

MANUFACTURED BLADES 





El 






MODE 

^NR> radians/sec 

cjr, radians/sec 

EXACT 

ESTIMATED 

EXACT 

ESTIMATED 

1 St 

2nd 

3rd 

17.3 

48.5 

95.5 

17.4 
50.0 
10 1.0 

49.2 
86.8 
1 37.0 

47.7 

85.7 
1 37.8 


1 st 

21.6 

21.1 

50.6 

49.2 

2nd 

58.9 

60.5 

92.4 

92.2 

3rd 

1 12.1 

1 22.0 

148.0 

154.0 


1 St 

21.9 

21.1 

74.0 

78.3 

2nd 

63.7 

59.5 

1 32.0 

1 34.4 

3rd 

126.0 

125.5 

200.0 

207.5 


1 st 

13.4 

14.6 

37.9 

37.8 

2nd 

43.7 

41.6 

71.0 

70.3 

3rd 

94.9 

94.5 

125.0 

124.0 


TIP ROOT TIP 



ROOT 
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TABLE II 


VALUES FOR INTEGRALS IN THE MODE-EXPANSION METHOD OF APPENDIX A 


p 


®PQo 


®pq.k 

®P<lt 

^<lo 





Hinged beams vith linear mass and stiffness distributions 


0 

1/3 

1/2 

1/4 

1 

0 

0 

1/3 

1/4 

r 


1 

0 

-0.18517 

-0.05578 

1 

0 

0 

0 

0.05577 

r 

0 

2 

0 

.09992 

-.00165 

1 

0 

0 

0 

.00165 

r 


3 

0 

-.06926 

-.00320 

1 

0 

0 

0 

.00321 

r 



0 

.05296 

-.00024 

1 

0 

0 

0 

.00025 

r 


5 

0 

-.04287 

-.00075 

1 

0 

0 

0 

.00077 

r 


1 

1.59938 

2.30532 

1.26052 

6.80791 

59.^3015 

25.63938 

1/4 

.14214 

r 


2 

- . 46528 

-1.22263 

-.18122 

3.59935 

0 

30.97991 

0 

.05286 

r 

1 

5 

-.09145 

.20921 

-.09117 

3.78925 

0 

- 4.62284 

0 

.00288 

r 


k 

-.03037 

-.21872 

-.02260 

3.85869 

0 

7.15092 

0 

.00531 

r 


5 

-.01293 

.12567 

-.01258 

3.8^21 

0 

-3.32571 

0 

.00082 

r 


2 

4.47610 

6.62225 

3 -43015 

17.79273 

624.12075 

299.60439 

1/4 

.12999 

r 

P 

5 

-1.54866 

-3.45753 

-.80926 

6.13170 

0 

239.3184 

0 

.05202 

r 


4 

-.57008 

.28431 

-.57093 

6.53295 

0 

-15.99519 

0 

.00182 

r 


5 

-.14453 

-.52531 

-.12984 

6.73955 

0 

48.65946 

0 

.00582 

r 


5 

8.99985 

-13.40607 

6.82541 

53-71970 

2716.90000 

1532.5410 

1/4 

.12741 

r 

5 

4 

-3.15890 

-6.67345 

-1.76094 

8.57761 

0 

892.5556 

0 

.05144 

r 


5 

-.80708 

.31838 

-.81443 

9.12863 

0 

-30.56181 

0 

.00150 

r 

ii 

4 

15.16836 

22.65800 

11.44593 

54.58150 

7945-0300 

3931.674 

1/4 

.12652 

r 


5 

-5.27717 

-10.91717 

-2.97035 

io. 9 &ro 7 

0 

2319.820 

0 

.05001 

r 

5 

5 

22.98186 

34.57920 

17.25665 

80.37810 

18500.2025 

9629.3965 

0 

.12710 

r 

Cantilever beams vith linear mass and stiffness distributions 


1 

0.29833 

0.39272 

0.23958 

1.16194 

3.09056 

0.59789 

1/4 

0.20163 

r 


2 

.17146 

.10558 

.18630 

1.84496 

0 

2.97335 

0 

.03838 

r 

1 

3 

-.19809 

-.26802 

-.13341 

.98538 

0 

-1.10249 

0 

.00508 

r 


4 

.13660 

.21828 

.09267 

1.64834 

0 

.93662 

0 

.00220 

r 


5 

-.11552 

-.19058 

-.07526 

1.14799 

0 

-.67480 

0 

.00094 

r 


2 

1.61955 

2.16178 

1.31905 

8.10433 

121.37958 

49.26172 

1/4 

.14854 

r 

p 

3 

-.04255 

-.47251 

.12583 

5.58811 

0 

64.85774 

0 

.04771 

r 

c. 

hr 

-.72797 

-.91085 

-.53287 

3.39561 

0 

-13.73693 

0 

.00514 

r 


5 

.47229 

.76571 

.30697 

5.70988 

0 

18.98574 

0 

.00430 

r 


3 

4.46488 

6.23803 

3.50050 

19.32474 

951.63772 

444.99392 

1/4 

.13308 

r 

3 

hr 

-.81857 

- 2.08457 

-.29328 

8.91208 

0 

567.5131 

0 

.04928 

r 


3 

-1.53861 

-1.78527 

-1.18157 

5.04050 

0 

-40.74667 

0 

.00317 

r 

k 

hr 

9.01587 

12.86480 

6.94974 

35.72554 

3654.3173 

1765.9083 

1/4 

.12905 

r 


3 

- 2.14253 

-4.75477 

-1.05264 

12.17992 

0 

1209.091 

0 

.05007 

r 

5 

3 

15.20032 

21.94810 

11.63995 

57-03549 

9985.9627 

4&70.7227 

1/4 

.12687 

r 
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TABLE III 


MODE RESULTS FOR NONROTATIKG CANTILEVER BEAMS WITH LINEAR 
MASS AND STIFFNESS DISTRIBUTIONS 


Station 

H H 
111 

7l' s 

•sl Hi 
III 

Y2 = 
($2 

Y2- s 
^2 

Y2" s 

9 S 2 " 

Hi 

III 

Hi 

III 

III 



m 

- 

El = EIo (exact solution, ref. 

20) 


0 

0 

0 

3.5160 

0 

0 

-22.0545 

0 

0 

61.6972 

1 

.0168 

.3274 

3.0352 

-.0926 

-1.6776 

-11.5406 

.2281 

5.7655 

14.0984 

2 

.0639 

.6065 

2.5508 

-.5011 

-2.3240 

-1.5432 

.6045 

3.1181 

24.3627 

3 

.1365 

.8378 

2.0775 

-.5261 

-2.0351 

6.9860 

.7562 

-.3551 

-40.5613 

k 

.2299 

1.0226 

1.6214 

-.6835 

- 1.0114 

12.9888 

•5259 

- 4.0599 

29.2300 

5 

.3395 

1.1631 

1.1958 

-.7137 

.4531 

15.7253 

.0197 

-5.5520 

1.2145 

6 

.4611 

1.2627 

.8085 

-.5895 

2.0194 

15.0599 

-.4738 

-3.7912 

32.4481 

7 

• 5959 

1.3266 

.4799 

-.3171 

3.3709 

11.5931 

-.657^ 

.3568 

46.6579 

1 8 

.7255 

1.3612 

.2246 

.0700 

4.2876 

6.6356 

-.5949 

4.735^ 

37.2963 

9 

.8624 

1.3745 

.0590 

.5258 

4.7095 

2 . 04 U 

.2285 

7.5585 

14.0713 

10 

1.0000 

1.3765 

0 

1.0000 

4.7808 

0 

1.0000 

7.8487 

0 




m = mo 

; El = EIo 




0 

1 

2 

3 

k 

5 

6 

7 

8 
9 

10 

0 

.0169 

.0642 

.1369 
.2304 
.3400 
.4617 
• 5914 

.7259 

.8626 

1.0000 

0.1695 

•4723 

.7271 

.9547 

1.0968 

1.2162 

1.2970 

1.5451 

1.3676 

1.3736 

I 

3.5104 

5.0282 

2.5482 

2.0760 

1.6207 

1.1937 

.8087 

.4806 

.2255 

.0595 

0 

0 

-.0925 

-.3008 

-.5257 

-.6830 

-.7131 

-,5890 

-.5166 

.0704 

.5239 

1.0000 

-0.9051 

-2.0829 

- 2.2495 

-1.5725 

-.3015 

1.2418 

2.7256 

3.8698 

4.5556 

4.7607 

-22.0247 

-11.5367 

-1.5444 

6.9810 

12.9814 

15.7169 

15.0516 

11.5857 

6.6278 

2.0377 

0 

0 

.2261 

.6024 

.7553 

.5262 

.0197 

-.4753 

-. 66 o 4 

-.3978 

.2270 

1.0000 

2.2608 

3.7628 

1.5296 

-2.2914 

-5.0647 

-4.9503 

-1.8512 

2.6262 

6.2482 

7.7299 

61.7316 

14.2680 

- 24.1735 

-40.5000 

-29.331^ 

1.0704 

32 . 4 i 48 

46.7572 

37.3956 

14.0565 

0 
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TABLE III.- Continiied 


MODE RJSULTS FOR NONROTATING CANTILEVER BEAMS WITH LINEAR 
MASS AND STIFFNESS DISTRIBUTIONS 


Station 

Yl 

Yl' 

Yi" 

Y2 

Y 2 ' 

^2 

Y5 

Y3' 

V " 

Y5 




m = 

mo; El 

: = Eio(i 

- 1 ) 




0 

1 

2 

5 

4 

5 

6 

7 

8 

9 

10 

0 

.0151 

.0582 

.1263 

.2159 

.525^ 

. 1(450 

.5769 

.7154 

.8572 

1.0000 

0.1508 

.4312 

.6807 

.8962 

1.0755 

1.2165 

1.5189 

1.5847 

1.4182 

1.4278 

3.0852 

2 . 8 o 44 

2.4942 

2.1551 

1.7910 

1.4102 

1.0259 

.6580 

.5550 

.0967 

0 

0 

-.0778 

-.2612 

-.4725 

-.6561 

-.6907 

-.5971 

-.5470 

.0555 

.5022 

1.0000 

-0.7779 

-1.8556 

-2.1111 

-1.6586 

-.5456 

.9562 

2.5002 

3.8252 

4.6675 

4.9778 

-17.9825 
-ID. 5955 

-2.7555 
4.8322 
11.1252 
15.0784 
15.9152 
13.4617 
8.4741 
2. 8758 
0 

0 

.l&Tl 

.5246 

.7053 

.5543 

.1097 

-5927 

-.6432 

-.4394 

.1846 

1.0000 

1.8715 

5.571+6 

1.8048 

-1.5085 

-4.4456 

- 5.0240 

-2.5050 

2.0377 

6.2401 

8.1542 

49.0088 
14.8158 
-16.8271 
-54.9664 
-31.1090 
-6.4431 
26.1221 
47.5240 
43.8272 
18. 6063 
0 



m 

= mo; El = EIo(l - 5t) 




0 

1 

2 

5 

k 

5 

6 

7 

8 
9 

10 

0 

.0125 

.0491 

.1088 

.1900 

.2908 

. 1(095 

.5426 

.6879 

.8417 

1.0000 

0.1247 

.5667 

.5966 

.8117 

1.0088 

1.1841 

1.5557 

1.4551 

1.5578 

1.5829 

2.5176 

2.4201 

2.2996 

2.1512 

1.9705 

1.7551 

1.4955 

1.1945 

.8469 

.4515 

0 

0 

-.0518 

-.1818 

-.5471 

-.4989 

-.5863 

-.5622 

-.3908 

-.0570 

.4244 

1.0000 

-0.5178 

-1.3001 

-1.6530 

-1.5179 

-.8741 

.2412 

1.7155 

5.5578 

4.8142 

5.7561 

-11.4951 

-7.8970 

-5.5774 

1.5516 

6.4680 

11.2476 

14.8955 

16.4929 

i 5 .o&rr 

9.7945 

0 

0 

.1079 

.5289 

.5018 

.4938 

.2551 

-. 157 *^ 

-.4768 

-4899 

.0240 

1.0000 

1.0789 

2.2098 

1.7288 

-.0799 

-2.3871 

-5.9249 

- 3 . 591(0 

-.1503 

5.1585 

9.7602 

26.4729 

11.4265 

-5.0568 

-18.7865 

-24.1585 

-16.5097 

4.6998 

55.2025 

54.9056 

49.9675 

0 
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TABLE III.- Continued 


MODE RESULTS FOR NONROTATING CANTILEVER BEAMS WITH LINEAR 
MASS AND STIFFNESS DISTRIBOTIONS 


Station 


Yl' 

Yi" 

^2 


V " 

^5 

^3' 

V 




m 

= 

- 1); El = EIo 




0 

0 


3.5961 

0 


-20.5149 

0 


56.8179 



0.1733 



-0.8526 



2.0307 


1 

.0173 


3.0782 

-.0853 


-10.2793 

.2031 


10.7772 



.4812 



-1. 8857 



3.1960 


2 

.0654 


2.5638 

-.2738 


-.6581 

.5227 


-24.7375 



.7373 



-1.9658 



.9179 


5 

.1592 


2.0614 

-.4704 


7.2497 

.6145 


-36.5657 



.9437 



-1.2637 



-2.5205 


k 

.2356 


1.5835 

-.5968 


12.4090 

.3624 


-22.1859 



1.1020 



-.0501 



-4.6174 


5 

.5438 


1.1442 

-.6018 


14.2921 

-0993 


6.8125 



1.2164 



1.5533 



-3.9663 


6 

.4654 


.7580 

-.4665 


13.0765 

-.4960 


32.2249 



1.2922 



2.6423 



-.8880 


7 

.5946 


.4391 

-.2022 


9.6161 

-.5848 


40.3627 



1.3362 



3.5963 



2.9901 


8 

.7283 


.1999 

.157^ 


5.2513 

-.2857 


29.5511 



1.3562 



4.1269 



5 - 8 r ^9 


9 

.8639 


.0509 

.5701 


1.5398 

.3018 


10.3404 



1.3613 



4.2990 



6.9825 


10 

1.0000 


0 

1.0000 


0 

1.0000 


0 




m = mo(^l - 1 ^; El = 


0 



0 

0 


3.1688 

0 


-16.8561 

0 


45.1241 



0.1546 



-0.7222 



1.6861 


1 

•0155 


2.8585 

-.0722 


- 9 . 53 l ^3 

.1686 


11.7757 



.4404 



-1.6728 



2.8950 


2 

.0595 


2.5168 

-.2395 


-1.8680 

.4581 


-17.8557 



.6921 



-1.8638 



1.2361 


5 

.1287 


2.1467 

-.4259 


5-2897 

.5817 


-32.5125 



.9068 



- 1.3480 



-1.8134 


k 

.2194 


1.7557 

-.5607 


10.8679 

.4004 


-24.9765 



1.0823 



-.2821 



-4.1661 


5 

.3276 


1.5565 

-.5889 


15.9565 

-.0162 


-.2658 



1.2180 



1.0866 



-4.1669 


6 

.4493 


.9647 

-.4802 


14.0175 

-4330 


27.5563 



1.3144 



2.4651 



-1.5214 


7 

.5808 


.6034 

-.2357 


11.3083 

-.5851 


42.1923 



1 . 37^+8 



3.5808 



2.5189 


8 

.7183 


.2979 

.1244 


6.7837 

-.3333 


55.4166 



1.4046 



4.2586 



5-9375 


9 

.8587 


.0836 

.5502 


2.1918 

-.2605 


13.9636 



1.4129 



4.4978 



7.3951 


10 

1 1.0000 


0 

1.0000 


0 

1.0000 


1.0000 
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TABLE III . - Continued 

MODE RESULTS FOR NONROTATING CAOTILEVER BEAMS WITE LINEAR 
MASS AND STIFFNESS DISTRIBUTIONS 


Station 


Yl’ 

Yi" 

^2 


Y ” 

^5 

^5’ 

Y 5 " 

“ = “o (l - 1 

); El = 

EIo(l - x) 



0 

0 


2.5984 

0 


- 11.0011 

0 


25.2707 



0.1284 



- 0.4911 



1.0099 


1 

.0128 


2.4797 

-.0491 


-7.2927 

.1010 


9-7964 



.3764 



- 1.2146 



1.9826 


2 

.0505 


2.5558 

-.1706 


-2.8931 

.2992 


-6.5750 



.6098 



-1.5002 



1.5592 


3 

.1115 


2.1566 

-.3206 


1.9575 

.4552 


-18.9042 



.8255 



- 1.3046 



-1.4506 


4 

.19to 


1.9455 

-.4511 


6.7820 

.5901 


-21.5881 



1.0200 



-.6320 



-2.5029 


5 

.2960 


1.6994 

-.5143 


10.9596 

.1598 


-11.5213 



1.1899 



.4505 



-5.5647 


6 

.4150 


1.4185 

-.4692 


15.7170 

-.2167 


9.5586 



1.5518 



1.8050 



- 2 . 6 o 4 i 


7 

• 5482 


1.1048 

-.2887 


14.4325 

-.4771 


35.1813 



1.4423 



5.2265 



.6360 


8 

.6924 


.7609 

.0339 


12.5408 

-.4155 


47.5995 



1.5184 



4.4562 



5.2112 


9 

.8445 


.5911 

.4796 


7.7223 

.1077 


59.7770 



1.5575 



5.2043 



8.9254 


10 

1.0000 


0 

1.0000 


0 

1.0000 


0 




m 

= mod 

- 5 f) ; El = EIo 




0 

0 


4.0558 

0 


-15.2162 

0 


35.5410 



0.1959 



-0.6069 



1.1740 


1 

.0194 


5.5450 

-.0607 


-6.2023 

.1174 


2.2191 



.5284 



- 1.2348 



1.4916 


2 

.0722 


2.6451 

-.1842 


1.8823 

.2666 


-19.6272 



.7929 



-1.0659 



-.2901 


3 

.1515 


1.9815 

-.2908 


7.6453 

.2575 


-19.7496 



.9911 



-.5292 



-2.1229 


k 

.2506 


1.3858 

-.5237 


10.2183 

.0255 


-2.7888 



1.1294 



.6672 



-2.3870 


5 

.3636 


.8803 

-.2570 


9.7558 

-.2155 


15.9480 



1.2175 



1.6238 



-.8854 


6 

.4853 


.4904 

-.0946 


7.2110 

-.3018 


23.7445 



1.2665 



2.3597 



1.5851 


7 

.6120 


.2218 

.1594 


4.0545 

-.1635 


18.8323 



1.2887 



2.7505 



3.2238 


8 

.T 4 o 8 


.0680 

.4144 


1.5122 

.1591 


8.5835 



1.2955 



2.9121 



4.1081 


9 

.8704 


.0068 

.7056 


.2282 

.5699 


1.4574 



1.2962 



2.9457 



4.3010 


10 

1.0000 


0 

1.0000 


0 

1.0000 


0 
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TABLE III.- Concluded 

MODE RESULTS FOR NONROTATING CANTILEVER BEAMS WITH LINEAR 
MASS AND STIFFNESS DISTRIBUTIONS 


Station 


Yl' 

V 

^2 

Y2' 

V " 

^5 

^5' 

Y " 

^5 




m = iiiq(1 - X 

); El = EIo(l - 1 

:) 



0 

0 


5.6228 

0 


- 12.6640 

0 


28,6145 



0.1752 



-0.5212 



0.9900 


1 

.0175 


5.1505 

-.0521 


-5.9369 

.0990 


3.6058 



.4903 



-1. n 48 



1.3986 


2 

.0666 


2.6356 

-.1636 


.8058 

.2389 


-15.6241 



.7559 



- 1.0442 



-.0314 


5 

.1419 


2.0961 

-.2680 


6.5475 

.2357 


-18.9598 



.9635 



-.4288 



-1.7917 


k 

.2383 


1.5602 

-.3109 


9.5678 

.0566 


- 6.0124 



1.1195 



.5049 



-2.3466 


5 

.5502 


1.0620 

-.2604 


10.0140 

-.1781 


12.4980 



1.2257 



1.4867 



-1.1598 


6 

.4728 


.6358 

-.1118 


8.1088 

-.2941 


25.4526 



1.2893 



2.2874 



1.0759 


7 

.6017 


.3106 

.1170 


4.9844 

-.1865 


21.2750 



1.3203 



2.7856 



5-1546 


8 

.7338 


.1054 

.3955 


2.0367 

.1270 


10.8683 



1-3307 



2.9997 



4.2345 


9 

.8668 


.0113 

.6955 


.3379 

.5504 


2.0489 



1.3318 



5.0448 



4.4958 


10 

1.0000 


0 

1.0000 


0 

1.0000 


0 




m = rao(l - x) 

) El = EIo(l - x) 



0 

0 


3.0750 

0 


-8.9965 

0 


18.1439 



0.1507 



-0.3854 



0.6710 


1 

.0151 


2.8299 

-.0385 


-4.9951 

.0671 


4.1930 



-4337 



-.8801 



1.0985 


2 

.0584 


2.5321 

-.1266 


-. 45&7 

.1769 


-8.8006 



.6869 



-.9247 



.2793 


5 

.0271 


2.1832 

-.2190 


4.0280 

.2049 


-14.5066 



.9052 



-5297 



-1.0798 


4 

.2176 


1.7937 

-.2720 


7.5654 

.0969 


-9.2434 



1.0846 



-.2128 



-1.9543 


5 

.3261 


1.3813 

-.2507 


9.4354 

-.0965 


4.3850 



1.2227 



1.1393 



-1.4924 


6 

.4484 


.9705 

-.1368 


9.2835 

-.2458 


18.4222 



1.5198 



2.0525 



.2825 


7 

.5803 


.5903 

.0684 


7.3013 

-.2175 


24.5923 



1.3788 



2.7735 



2.6264 


8 

.7182 


.2735 

•3458 


4.2407 

.0451 


18.9564 



1.4061 



3 . 19 &T 



4.4100 


9 

.8589 


.0547 

.6657 


1.3128 

.4861 


6.8571 



1 . 4 n 6 



5.3434 



5.1389 


10 

1.0000 


1.0000 

1.0000 


0 

1.0000 


0 


J 
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TABLE IV 

MODE RESULTS FOR NONROTATING HINGED BEAMS WITH LINEAR MASS AND STIFFNESS DISTRIBUTIONS 


Station 

H 

Ill 

”h 

Ill 

III 

III 

III 

”c\J 

Ill 

III 

Yj' s (ilj' 

Yj" s 




m = mo; 

El = EIo 

(exact solution, ref. 20) 



0 

1 

0 

-2.7002 

0 

0 

5.0043 

0 

0 

-7.2193 

0 

2 

3 

k 

-,4830 

-1.8617 

7.9756 

.7001 

.7950 

-34.8133 

-.6299 

3.2792 

65.6943 

5 

6 
7 

-.6620 

.1938 

11.6061 

.2257 

-4.7026 

-10.5596 

.5732 

4.2548 

-59.5223 

8 

9 

-0973 

2.3756 

9.3030 

-.6005 

-2.0600 

32.9576 

.1190 

- 7.0448 

-10.6536 

10 










11 










12 

.2274 

3.6749 

3.5134 

-.2940 

4.9033 

26.8434 

-.6076 

2.8933 

76 .&T 02 

13 










14 










15 

1.0000 

3.9297 

0 

1.0000 

7.0686 

0 

1.0000 

10.2102 

0 





m = mo; 

El = EIo 





0 

1 

2 

3 

4 

5 

6 
T 
8 
9 

10 

11 

12 

13 

14 

15 

0 

•1778 

-.3428 

-.4830 

- 5 &T 9 

-.6494 

-.6620 

-.6233 

-.5339 

-0973 

-.2195 

-.0083 

.2274 

.4788 

.7382 

1.0000 

-2.6675 

-2.4751 

-2.1026 

-10737 

-.9222 

-.1890 

0809 

1.3413 

2.0490 

2.6669 

3.1679 

3.5357 

3.7713 

3.8904 

3.9269 

0 

2.9017 

5.6183 

7.97^8 

9.8240 

11.0559 

11.6059 

11.4610 

10.6619 

9.3030 

70279 

5.5239 

3.5133 

1.7444 

.4821 

0 

0 

.3217 

.5731 

.6997 

. 67 i ^5 

.5034 

.2246 

-.1000 

-.3977 

-.6006 

-.6595 

-.5667 

-.2935 

.0847 

.5309 

1.0000 

4.8253 

3.7708 

1.9000 

-.3784 

-2.5660 

-4.1825 

-4.8692 

- 4.4659 

-3.0451 

-.8832 

1 . 5 &T 8 

5.9028 

5.6720 

6.6930 

7.0367 

0 

-16.1228 

-28.5821 

-54.8100 

-33.^256 

-24.7037 

-10.5043 

6.1398 

21.70^ 

32.9430 

57 - 659 V 

55-2256 

26.7997 

15.2316 

if .6659 

0 

0 

-.4342 

-.6814 

-.6298 

-.3017 

.1561 

.5389 

.6758 

.5070 

.1106 

-.3311 

-.6112 

-.5878 

-.2400 

.3365 

1.0000 

-6.5157 

-3.7078 

. 77‘^7 

4.9212 

6.8673 

5.7417 

2.0541 

-2.5328 

- 5.9464 

-6.6245 

- 4.2021 

.3507 

5.2169 

8.6477 

9.9528 

0 

43.4748 

69.7829 

64.4864 

30.1462 

-17.7598 

-57.7175 

-71.7578 

-53.4956 

-10.8931 

37.3827 

70.4947 

75.1702 

52.3442 

18.2343 

0 
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TABLE IV.- Continued 

MODE RESULTS FOR NONROTATING HINGED BEAMS WITH LINEAR MASS AND STIFFNESS DISTRIBUTIONS 


Station 



Yl" 

T^2 

■^2 


^5 

^5’ 

V 





m = eoq; 

El = EIo^l 

2 1 




0 

1 

2 

3 

k 

3 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

0 

-.1601 

-.3104 

-.4411 

-.5430 

-.6079 

-.6293 

-.6030 

-.5272 

-.4031 

-.2347 

-.0283 

.2076 

. 4^35 

.7300 

1.0000 

-2.1(011 

-2.25U5 

-I.96ID 

-1.5283 

-.9752 

-.3211 

.3949 

1.1363 

1.8610 

2.5269 

3.0959 

3.5582 

5-8577 

3.9981 

4.0501 

0 

2.1968 

4.4128 

6.5121 

8.3585 

9.8221 

10.7894 

11.1731 

10.9231 

10.0375 

8.5728 

6.6547 

4.4883 

2.3693 

.6982 

0 

0 

.2825 

.5112 

.6413 

.6441 

.5151 

.2771 

-.0227 

-.5201 

-.5466 

-.6433 

-.5747 

-.5369 

.0406 

.50^7 

1.0000 

4.2349 

3.4337 

1.9512 

. 04 l 4 

-1.9542 

-3.5700 

-4.4964 

- 4.4618 

- 3-3977 

-1.4512 

1.0303 

5.5665 

5.6618 

6.9620 

7.4296 

0 

-12.1633 - 

-22.5509 

-29.0968 

-30.1375 

-24.9921 

- 14.2022 

.4449 

16.1636 

29.6539 

37-8479 

38 . 6 &T 4 

31.8978 

19.5720 

6.4608 

0 

-.3806 

-.6187 

-.6118 

-.3523 

.0612 

.4585 

.6662 

.5843 

.2385 

-.2199 

-.5723 

-.6286 

-.3172 

.2820 

1.0000 

-5.7084 

-5.5718 

.0330 

3.8918 

6.2024 

5.9598 

3.1149 

-1.2129 

-5.1621 

-6.8767 

-5.2857 

-.8435 

4.6710 

8.9872 

10.7702 

0 

52. 81(85 
57-1675 
59.0998 

36 .J 01 l( 

- 3 - 51 ( 1(5 

•45.9333 

-67.5091 

-61.8610 

-26.7887 

a 4.2333 

68.7551* 

85.5598 

66 . 4741 ( 

25.3995 

0 



m = BOo; 

El = EIo(l 

- S) 




0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

n 

12 

13 

14 

15 

0 

-.1241 

-.2425 

-.3492 

-.4377 

-.5015 

-.5344 

-.5308 

-.4860 

-.3966 

-.2609 

-.0795 

.1447 

.4059 

.6951 

1.0000 

-1.8612 

-1.7765 

-1.6000 

-1.3272 

-.9572 

-.4940 

.0539 

.6722 

1.3411 

2.0346 

2.7212 

3.5654 

3.9182 

4.3585 

4.5750 

0 

1.2599 

2.6412 

4.0919 

5.5538 

6.9600 

8.2356 

9.2994 

10.0650 

op. 4435 
10.5457 
9.6854 
8.3825 
6.3673 
5.5845 
0 

0 

.1893 

.3518 

.4622 

.4997 

.4526 

.3215 

.1222 

-.1135 

-.3^2 

-.5027 

-.5451 

-.4204 

-.1047 

.3896 

1.0000 

2.8391 

2.4377 

1.6557 

.5632 

-.7071 

-1.9666 

-2.9887 

-3.5362 

-3.5992 

-2.4389 

-.6351 

1.8695 

4.7359 

7.4i4o 

9.1568 

. 

0 

-6.01403 

-11.8108 

-16.5435 

-19.2804 

-19.1626 

-15.6121 

-8.4580 

1.9007 

14.3995 

27.2556 

58.0134 

45.7041 

41.1267 

27.2670 

0 

0 

-.2299 

-.3984 

-.4489 

-.5552 

-.1535 

.1449 

.3882 

.4954 

.5964 

.1086 

-.2615 

-5209 

-.4517 

.0822 

1.0000 

- 3.4482 

-2.5277 

-.7583 

1.4065 

3.2949 

4.2065 

3.6496 

1.5779 

-1.4555 

-4.3171 

-5.5508 

-5.8912 

1.0372 

8.0094 

13.7667 

0 

15.9649 

27.3446 

55.6311 

29.5841 

14.6758 

-7.9710 

-31.5748 

- 46.9233 

- 44.9846 

-20. 5659 
23.7570 
74.7164 
108.0221 
92.3208 
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TABI£ IV.- Continued 

MODE RESULTS FOR NONROTATING HINGED BEAMS WITH LINEAR MASS AND STIFFTTESS DISTRIBUTIONS 


Station 



V 

^2 

^^2 

V 

^5 

^3' 

^3” 





m = ihq ^1 

- 1 ); El = : 

EIo 




0 

1 

2 

3 

k 

5 

6 

T 

8 

9 

10 

11 

13 

14 

15 

0 

-.1558 

-.2990 

-.4182 

-.5038 

-.5^^84 

-.5478 

-.5005 

-.4083 

-.2751 

-.1071 

.0884 

.3037 

.5312 

.7647 

1,0000 

-2.5363 

-2.1488 

-1.7886 

-1.2830 

-.6693 

.0096 

.7083 

1.3836 

1.9975 

2.5203 

2.9330 

3.2285 

5.4124 

3.5029 

3.5298 

0 

2.8306 

5.^383 

7.6302 

9.2606 

10.2394 

10.5356 

10.1775 

9.2456 

7.8665 

6.1992 

4.4254 

2.7358 

1.3192 

.3538 

0 

0 

.2727 

.4785 

.5674 

.5186 

.3447 

.0&T2 

-.1932 

-.4314 

-.5716 

-.5787 

-.4431 

-.1812 

.1730 

.5782 

1.0000 

4.0912 

3.0857 

1.3535 

-.7308 

-2.6093 

-3.8621 

-4.2068 

-5.5727 

-2.1027 

-.1044 

2.0318 

3.9288 

5.3124 

6.0784 

6.3271 

0 

-15.4147 

-26.8539 

-31.5958 

-28.7523 

-19.1^^53 

-5.2516 

9.6994 

22.4495 

30.4780 

32.5262 

28.8072 

20.8930 

11.3644 

3.5416 

0 

0 

-.3730 

-.5607 

-.4680 

-.1402 

.2654 

.5497 

.5846 

.3524 

-.0405 

-.4169 

-.6068 

-.5150 

-.1530 

.3921 

1.0000 

-5.5947 

-2. 6156 

1.3899 

4.9171 

6.0537 

4.2957 

.5225 

-3.4830 

-5.8954 

-5.6463 

-2. 8473 

1.3765 

5.4594 

8.1471 

9.1186 

0 

43.3868 

65.7475 

55.0307 

17.6072 

-27.57^5 

-58.9163 

-62.3944 

-37.4169 

3.9909 

43.5972 

65.5174 

62.9664 

40.7963 

13.4256 

0 




m : 

IXICM 

1 

');£!= Elofl - 1 ) 




0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

0 

-.1405 

-.2714 

-.3831 

-.4670 

-.5157 

-.5236 

-.4&rr 

-.1407^^ 

-.2847 

-.1240 

.0684 

.2847 

.5168 

.7571 

1.0000 

-2.K>75 

-1.9632 

-1.676U 

-1.2583 

-.7296 

-.1191 

.5384 

l.ao46 

1.8407 

2.4110 

2.8856 

5.2447 

5.4812 

5.6044 

5.6433 

0 

2.1650 

4.3162 

6.2980 

7.9670 

9.2014 

9.9102 

10.0412 

9.5&T5 

8.5912 

7.1455 

5.3954 

3.5362 

1.8125 

.5182 

0 

0 

.2937 

.4281 

.5231 

.5013 

.3642 

.1395 

-.1251 

-.369^ 

-.551^7 

-.5760 

-.4711 

-.2251 

.1318 

.5545 

1.0000 

3.5959 

2.8262 

1.4247 

-.3274 

-2.0558 

-3.3735 

-3.9664 

-3.6638 

-2.4809 

-.6194 

1.57^ 

5.6896 

5.3528 

6.5^107 

6.6828 

0 

-11.7159 

-21.3762 

-26.7516 

-26.4302 

-20.15&7 

-9.1019 

4.5670 

18.0105 

28.3717 

33.4411 

32.2076 

25.2411 

14.7929 

4.6798 

0 

0 

-.3245 

-.5079 

-.4640 

-.2030 

.1630 

.4700 

.57^ 

.4204 

.0681 

-.3262 

-.5765 

-.5468 

-.2092 

.3517 

1.0000 

-4.8668 

-2.7515 

.6587 

3.9139 

5.*^909 

4.6057 

1.5685 

-2.5134 

-5.2845 

-5.91‘^9 

-3.75*^5 

.4454 

5.0645 

8.4144 

9.7239 

0 

32.6840 

53.5102 

50.6902 

24.6040 

-15.7307 
-47. 2762 
-60.4825 
-46.3612 
-9.9582 
33.4867 

65.2075 

71.5806 

51.3484 

18.4557 

0 
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TABLE IV.- Continued 

MODE RESULTS FOR NONROTATING BINGED BEAMS WITH LINEAR MASS AND STIFINESS DISTRIBUTIONS 


Station 

Yl 

TCl' 

V 

^2 

^2' 

V 

^5 

^ 3 ' 





m 

= “0(1 - 1 

.); El = Elod - 5!) 

0 

0 


0 

0 


0 

0 


0 



-1.6543 



3.4651 



-3.0051 


1 

-.1103 


1.2820 

.1643 


-6.0465 

-.2003 


14.3350 



-1.5682 



2.0642 



-2.0645 


2 

-,2XkQ 


2.6679 

.3020 


-11.6724 

-.3380 


26.9253 



-1.3902 



1.2932 



-.3264 


5 

-.5075 


4.0886 

.3882 


-15.9371 

-.3597 


31.1310 



-1.1178 



.2431 



1.6735 


k 

-.5500 


5.4721 

.i4044 


-17.8592 

-.2482 


24.3216 



-.7556 



-.9314 



3.2227 


5 

-.4523 


6.7435 

.5423 


-16.7340 

-.0333 


7.4502 



-.3051 



-2.0293 



3.6744 


6 

-. 1 i -526 


7.8273 

.2070 


-12.2504 

.2116 


-14.6474 



.2153 



-2.8297 



2.6991 


7 

-A583 


8.6503 

.0184 


-4.6272 

.3916 


-34.2032 



.7902 



-3.1264 



.4730 


8 

-.3856 


9.1450 

-.1901 


5.3591 

.4231 


-42.7143 



1.3977 



-2.7646 



-2.2779 


9 

-.292lf 


9.2516 

-.3744 


16.3886 

.2712 


-33.8828 



2.0120 



-1.6775 



-4.4252 


10 

-.1583 


8.9210 

-.4862 


26.6827 

-.0238 


-6.5779 



2.6o4i 



.0846 



-4.7776 


n 

.0153 


8.1173 

-.4806 


54.1999 

-.3423 


33.1451 



3.1425 



2.3363 



-2.5501 


12 

.2248 


6.8194 

-.3248 


56.8503 

-.5123 


71. 8610 



3.5944 



4.7546 



2.1574 


13 

.4645 


5.0225 

-.0078 


52.7806 

-.3684 


90.9^56 



3.9265 



6.8945 



8.0186 


Ik 

.7262 


2.7589 

.4518 


20.6465 

.1667 


71.1259 



4.1066 



8.2252 



12.4991 


15 

1.0000 


0 

1.0000 


0 

1.0000 


0 





m = mo(l 

- x); El = 

EIo 




0 

0 


0 

0 


0 

0 


0 



-1.5091 



2.1411 



-2.5582 


1 

-.1006 


2.6772 

.1427 


-11.3779 

-.1705 


27.6738 



-1.3324 



1.4053 



-.8148 


2 

-.1894 


5.0299 

.2364 


-18.6533 

-.2249 


36.4609 



-1.0003 



.1972 



CVJ 

H 


3 

-.2561 


6.8123 

.2496 


-19.5232 

-.1258 


20.8593 



-.5501 



-1.0696 



2.8127 


k 

-.2928 


7.8813 

.1783 


-14.1242 

.0617 


-7.4162 



-.0288 



-1.9884 



2.35^9 


5 

-.2947 


8.2013 

.0457 


-4.5930 

.2isr 


-30.0183 



.5141 



-2.2893 



.4591 


6 

-.2604 


7.8349 

-.1069 


5.95^8 

.2493 


-34.1493 



1.055^ 



-1.9048 



-1.7115 


7 

-.1916 


6.9215 

-.2359 


14.4825 

.1352 


-19.1967 



1.4928 



-.9611 



-2.9^17 


8 

.0920 


5.6493 

-.2980 


19.1285 

-.0409 


5.1045 



1.8686 



.2897 



-2.6244 


9 

.0525 


4.2248 

-.2787 


19.3997 

-.2358 


25.9492 



2.1505 



1.5654 



-.9652 


10 

.1759 


2. 8430 

-.1744 


16.1429 

-.3002 


34.8202 



2. 3411- 



2.6292 



1.2810 


11 

.5320 


1.6634 

.0008 


11.0251 

-.2148 


30.8241 



2.4537 



3.3640 



3.2917 


12 

.4956 


.7904 

.2251 


5.&T86 

.0047 


19.3206 



ro 

$5 



3.7635 



k.^nQ 


15 

.6628 


.2606 

.4760 


2.1146 

.3098 


7.7297 



2.5274 



3.9151^ 



5.1208 


Ik 

.8315 


.0559 

.7370 


.3109 

.6512 


1.2214 



2.5308 



3.9444 



5.2315 


15 

1.0000 


0 

0 


0 

1.0000 


0 
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TABLE rV.- Concluded 


MODE RESUUS FOR NONROTATING HINGEID BEAMS WITH LINEAR MASS AND STIFFNESS DISTRIBUTIONS 


Station 

Yl 

II' 

Yl" 

Y 2 

Y2' 

Y2" 

^3 

^3' 

V 




m = 

®o(l - x); 

El = EIo ^1 

-I) 



0 

1 

2 

3 

4 

5 

6 
T 
8 
9 

10 

11 

12 

13 

14 

15 

0 

-.0912 

-.1731 

-.2366 

-.2743 

-.2808 

-.2535 

-.1922 

-.0995 

.0203 

.1615 

.3179 

.4838 

.6545 

.8271 

1.0000 

-1.3685 

-1.2273 

-.9526 

-.5655 

-.0984 

.4100 

.9190 

1.3913 

1.7972 

2.UL76 

2.3460 

2.4884 

2.5613 

2.5883 

2.5934 

0 

2.1268 

4.1481 

5.8470 

7-0565 

7-6765 

7.6819 

7.1221 

6.1118 

4.8141 

3.4186 

2.1148 

1.0648 

.5729 

.0547 

0 

0 

.1254 

.2125 

.2332 

.1796 

.0658 

-.0759 

-.2045 

-.2813 

-.2802 

-.1924 

-.0264 

.1975 

.4550 

.7261 

1.0000 

1.8817 

1.3053 

.3843 

-.8044 

-1.7075 

-2.1257 

-1.9280 

-1.1527 

.0165 

1.3170 

2.4906 

5.3583 

5.8625 

4.0667 

4.1079 

0 

- 8.8479 

-15.2839 

-17.1710 

-15.8954 

- 6.4412 

3.0165 

11.8656 

17 .&T 50 

19.8367 

17-8332 

13.0866 

7.4790 

2.8817 

.4541 

0 

0 

-.1505 

-.2086 

-.1360 

.0265 

.1836 

.2419 

.1637 

-.0154 

-.1976 

-.2913 

-.2359 

-.0314 

.2773 

.6335 

1.0000 

-2.2572 

-.8715 

1.0885 

2.4377 

2.3568 

.&T 50 

-1.1742 

-2.6558 

-2.7629 

-1.4052 

.8298 

5.0677 

4.6312 

5.3430 

5.4972 

0 

21.7677 

30.8281 

21.1573 

-1.3903 

-25.5985 

- 32.2246 

-23.1981 

-1.5876 

21.2899 

54.8049 

34.5490 

23.7284 

10.2955 

1.7510 

0 




m = 

= ni^(l - X 

); El = EIo(l - x) 




0 

1 

2 

3 

4 

5 

6 
T 
8 
9 

10 

U 

12 

13 

14 

15 

0 

-.0748 

-.1433 

-.1989 

-.2^56 

-.2480 

-.2323 

-.1862 

-.1095 

-.0036 

.1282 

.2813 

.4502 

.6275 

.8139 

1.0000 

-1.1222 

-1.0267 

-i.25&r 

-.5499 

-.1864 

.2356 

.6908 

1.1511 

1.5883 

1.9768 

2.2964 

2.5545 

2.6884 

2.7669 

2.7912 

0 

1.4218 

2.8943 

4.2872 

5.4778 

6.3621 

6.8635 

6.9403 

6.5903 

5.8551 

4.8088 

3.5747 

2.2986 

1.1517 

.32)6 

0 

0 

.0915 

.1602 

.i&rd 

,1620 

. 0 &T 9 

-.0193 

-.1331 

-.2216 

-.2553 

-.2124 

-.0846 

.1218 

.3878 

.6882 

1.0000 

1.3728 

1.0305 

.4020 

-.3755 

-1.1116 

-1.6078 

-1.7061 

-1.3297 

-.5045 

.6438 

1.9177 

3.0957 

3.9904 

4.5050 

4.6775 

0 

-5.2011 

-9.6016 

-11.9049 

-11.2941 

- 7.6441 

-1.5725 

5.6793 

12.5389 

17.4785 

19.3913 

17.9102 

13.5472 

7.6745 

2.5573 

0 

0 

-.1037 

-.1562 

-.1277 

-.0274 

.0975 

.1824 

.1774 

.0741 

-.0852 

-.2212 

-.2503 

-.1201 

.1683 

.5645 

1.0000 

-1.5560 

-7&n 

.4284 

1.5050 

1 .&T 22 

1.2737 

-.0748 

-1.5500 

-2.3886 

-2.0398 

-.4364 

1.9527 

4.3263 

5.9419 

6.5395 

0 

11.8732 

18.9054 

16.7673 

5.7258 

-9.3345 

-21.0490 

-23.0546 

-15.1568 

3.3345 

24.8731 

37.0753 

36.7275 

24.6765 

8.3232 

0 


TABLE V 


VJl 

00 


MODE (XEFFICIENTS FOR ROTATING BEAMS HINGED AT THE ROOT 


5 

yn = 

q=0 




Anq 








Linear mass and 




TMiform mass and stiffness distributions 


stiffness distributions 

n 

q 






(®t = = 0 ) 



H 

II 




“/“TO; 

L=1 



5 = 0 ^ 

5 = 10 ^ 

0 

II 

Id) 

S = 0 ^ 

e = 10 ^ 



r = 

0 

r = 0.1 

r = 0.5 

r = 1 

r = 

= 0 


0 

0 

-0.00676 

-2.99989 

-1.50000 

-3.00000 

0.38847 

0.38042 


1 

.90279 

.88995 

3.59557 

2.19920 

3.57296 

.57936 

.57921 


2 

.08290 

. 09^82 

.30674 

.22694 

.32128 

.02942 

.03739 

1 

3 

.01170 

.01223 

.06637 

.05056 

.07153 

.00226 

.00234 


k 

.00209 

.00150 

.02200 

.01656 

.02408 

.00053 

.00071 


5 

.00052 

.00827 

.00920 

.00673 

.01011 

-.00003 

-.00007 


0 

0 

.00111 

-3.00008 

-1.50000 

-3.00000 

.21916 

.21681 


1 

-.08920 

-.10460 

-5.12037 

-2.72581 

-5.17981 

.31280 

.30752 


2 

.95499 

.95733 

7.59936 

4.20150 

7.51089 

.43608 

.43748 

2 

3 

.11344 

.13274 

1.08367 

.75215 

1.18259 

.02899 

.05480 


k 

.01745 

.02015 

. 312 &T 

.21113 

.34805 

.00216 

.00241 


5 

.00352 

-.00672 

.12455 

.08096 

.13835 

.00082 

.00098 


0 

0 

-.00028 

-3.00002 

-1.50000 

-5.00000 

.15969 

.15854 


1 

-.00209 

.00029 

-4.11650 

- 2.04984 

-4.10909 

.21031 

.20886 


2 

-.12027 

-.14548 

-6.70530 

-3.73548 

-6.98322 

.23822 

.25474 


3 

.99297 

1.02956 

12.14772 

6.54922 

12.06709 

.36653 

.36824 

5 

4 

.11228 

.13716 

2.01361 

1.51803 

2.27887 

.02448 

.02874 


5 

.01711 

-.02126 

.66071 

. 4 l 8 o 4 

.74647 

.00077 

.00089 
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TABLE VI 


MODE COEFFICIENTS FOR ROTATING CANTILEVER BEAMS 


5 

yn = ZI 

q=l 




N 

n 



Ha 

Iform mass and stiffness 

distributions 




Linear mass and 
stiffness distributions 
(int = Elt = 0 ) 



“/“to; 

= 2 
L 

“/“to; 


^/“to 


^/<%R 

1=^ 

“/“toi “ 

7.37 

“/“tox ■ 

10.43 

“/“tox 

10.42 

“/“tox " 

14.76 

“/“tox 

- 1.97 

“/“tox 

= 5.91 



5 = 0 ^ 

e = 10 ^ 

e = 0 ^ 

e - 10 ^ 

e = 0^ 

e = 10 ^ 

e = 0 ^ 

e = 10 ^ 


e = 



e = 0 ^ 

e = 10 ^ 

e = 0 ^ 

e = 10^ 





r = 

= 0 

r = 0.5 

r = 1 

r = 0.5 

r = 1 


r = 

= 0 


1 

1 

2 

3 

4 

5 

I.OI4O97 

-Olt -833 

.00871 

-.00204 

.00069 

1.04089 

-.04894 

.00955 

-.00230 

.00080 

1.06395 

-.07740 

.01624 

-.00445 

.00165 

1.06197 

-.07629 

.01732 

-.00486 

.OOl&T 

1.07930 

-.09825 

.02318 

-.00713 

.00290 

1.07532 

-.09511 

.02415 

-.00756 

.00320 

1.09529 

-.12241 

.03332 

-.01177 

.00557 

1.08853 

-.11612 
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Figure 1.- Beams treated by both the "exact” and Rayleigh methods 
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Figure 2.- Effect of rotational speed on the bending frequencies of a 

■uniform hinged beam. 
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Figure 5.- Effect of rotational speed on the bending frequencies of a 
hinged beam with linear mass and stiffness distribution. 
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Figure k. - Con5)arison of frequencies of uniform and "linear" Mnged beams 

with zero offset. 
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Figure 5 .- Effect of rotational speed on the bending frequencies of a 

uniform cantilever beam. 
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Figure 6 .- Effect of rotational speed on the bending frequencies of a 
cantilever beam with linear mass and stiffness distribution. 
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Figure 7.- Comparison of frequencies of uniform and "linear" cantilever 

beams with zero offset. 
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Figure 



8. - Bending frequencies of a uniform cantilever beam with a mass 
at the tip and with zero offset, r = 

Beam mass 
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Figure 9.- Bending frequencies of a uniform hinged beam with a mass at 

Tip mass 

the tip and with zero offset, r = — . 

Beam mass 
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Figiare 10. - Bending frequencies of a \xnifonn hinged beam with a mass at 
the tip and with zero offset as a f\mction of cantilever -beam 

rotational-speed parameter, r = — mass ^ 

Beam mass 
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Figure 11.- Bending frequency coefficients for hinged beams with 

linear mass and stiffness distributions. 
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Figure 12.- Zero-offset Southwell coefficient KA for hinged beams 
with linear mass and stiffness distributions. 
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Figure IJ.- Offset-correction factors for 
for hinged beams with linear mass and 


Southwell coefficients 
stiffness distributions 
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Figure li<-. - 


Bending frequency coefficients an for cantilever beams 
with linear mass and stiffness distributions. 


L 



NACA TN 5^^-59 


7I+ 





Figure 15 .- Zero-offset Southwell coefficients Kq^ for cantilever beams 
with linear mass and stiffness distributions- 
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Figure l6.- Offset-correction factors for Southwell coefficients 

for cantilever beams with linear mass and stiffness distributions. 
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Figure 17* - Bending frequency coefficients for nonrotating uniform canti- 
lever beams with a mass at the tip. r = 


Beam mass 
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Figure l8.- 


Zero-offset Southwell coefficients for a uniform cantilever 
beam with a mass at the tip. r = mas ^^ 


Beam mass 
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Figure 19.- Bending frequency coefficients for nonrotating \miform canti- 
lever beams with a mass at the tip. r = ^ 

Beam mass 
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Figure 20.- 


Zero-offset Southwell coefficients for a uniform hinged beam 


with a mass at the tip. r 


Tip mass 
Beam mass 
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Figure 21.- Offset-correction factors for Southwell coefficients for the 
pendulum mode of hinged beams with linear mass distribution plus a 
mass at the tip. 
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(a) Cantilever beams with linear stiffness distribution. 

Figure 22 .- A new rotating beam frequency coefficient which is essentially 
independent of beam mass distribution. 
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(b) Hinged beams with linear stiffness distribution. 
Figure 22. - Concluded. 
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Figure 23.- Bending inodes of nonrotating hinged beams with linear mass 

and stiffness distributions. 
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Figure 24.- Bending inodes of nonrotating cantilever beams with linear 

mass and stiffness distributions. 
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Figure 25 . - Comparison of bending modes of a rotating and nonrotating 

uniform hinged beam. 
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Figure 26.- Comparison of bending modes for a rotating and nonrotating 
hinged beam vith linear mass and stiffness distribution. 
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Figure 27-- Coniparison of bending modes of a rotating and nonrotating 

uniform cantilever beam. 
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Figirre 28.- Comparison of bending modes of a rotating and nonrotating 
cantilever beam with linear mass and stiffness distribution. 
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"Figure 29-- Comparison of bending modes of a rotating and nonrotating 
xmiform hinged beam with a tip mass equal to the mass of the beam. 
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Figure 50- - Comparison of bending inodes of a rotating and nonrotating 
uniform cantilever beam with a mass at the tip. 
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